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ABSTRACT 


Two  programs  are  described  which  have  been  coded  for  the 
IBM-704  (and  are  compatible  with  the  IBM-7090).  The  programs 
calculate  the  neutron  and/pr  gamma  spectra,  heat  generation  rate, 
and/or  dose  rate  at  each  of  a  group  of  point  detectors,  due  to 
each  of  a  group  of  point  sources.  The  sources  may  be  divided  Into 
sets,  with  each  set  having  a  unique  source  spectra.  In  addition 
to  the  above  calculation,  the  spectrum,  heating  rate,  and/or 
dose  rate  for  each  detector,  summed  over  each  source-point  set 
and  over  the  entire  source  group  may  be  computed.  The  two  pro¬ 
grams  are  similar,  both  computationally,  and  as  regards  Input  and 
output  Information,  excepting  for  the  complexity  of  the  problem 
geometries  acceptable  to  the  programs. 
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I  INTRODUCTION 


The  two  shield-penetration  programs  (C-17  and  L-63)  described 
here  were  coded  to  take  advantage  of  the  data  given  In  Reference  1 
as  well  as  to  allow  the  computation  of  the  heat  generation  rates  In 
a  shield  system  due  to  the  fast-neutron  and  gamma-ray  reactor  leak¬ 
age  spectra. 

These  programs  are  the  latest  generation  In  the  continued 
development  of  moments-method  shield-penetration  programs  conducted 
at  the  Nuclear  Aerospace  Research  Facility  (NARP)  at  General 
Dynamlcs/Fort  Worth  (GD/FW).  These  methods  were  first  used  at 
NARF  to  analyze  data  from  the  Nuclear  Test  Aircraft  flights.  The 
methods  were  also  used  to  code  shield-penetration  programs  for  the 
JBM-701  and  IBM-704  In  order  to  compute  fast-neutron  and  gamma-ray 
spectra  and  dose  rates  In  shield  systems  described  by  a  geometry 
system  similar  to  the  simple  geometry  routine  described  In  Section  II. 
This  family  of  programs  has  been  developed  around  the  following 
basic  concepts: 

Moments  Method  -  This  method  Is  based  upon  the  differential 
energy  spectra  calculated  by  the  Nuclear  Development  Corporation 
of  America  (NDA)  for  a  point  Isotropic  source  In  an  Infinite 
medium  (Ref.  2).  The  data  represents  a  moments-method  solution 
of  the  fast-neutron  and  gamma-ray  transport  equation.  A  detailed 
description  of  the  use  of  these  data  Is  given  In  Reference  3« 

This  method  Implies  that  the  only  portion  of  the  system  affecting 
dose  rate  (or  spectra,  or  heat  generation  rate)  at  a  detector  due 


11 


to  a  point  source  a  distance  r  from  the  detector  Is  that  portion 
of  the  system  lying  on  the  line  of  sight  between  the  source  and 
the  detectorj  hence,  the  techniques  discussed  below  are  required 
In  the  use  of  this  method. 

The  Stepping  Point  Method  -  This  Is  an  Iterative  method  of 
determining  the  Intercepted  distances  between  two  points.  In  this 
procedure,  an  Iterative  scheme  Is  used  to  "step"  a  point  through 
the  volumes  of  the  system,  and,  after  each  step,  tests  are  performed 
to  determine  whether  a  boundary  of  the  volume  has  been  crossed. 

The  Intercepted  distance  Is  then  determined  by  the  number  of  steps 
and  the  length  of  each  step  taken  In  each  volume.  A  discussion 
of  the  basic  Ideas  and  methods  used  In  adapting  this  concept  to 
computer  programs  Is  given  In  Reference  4. 

Distributed  Source  -  The  moments-method  data  used  In  the 
calculations  are  for  point  Isotropic  sources;  thus.  It  Is 
necessary  to  approximate  the  leakage  from  a  reactor  or  other 
source  by  a  set  of  point  sources.  Each  of  these  sources  repre¬ 
sents  the  radiation  born  In  an  elemental  volume  containing  a 
point  so  that  an  Integration  of  source  points  over  the  source 
volume  must  equal  the  total  radiation  generated  by  the  source. 

The  new  programs  differ  from  the  earlier  versions  princi¬ 
pally  In  the  following  new  features? 

1.  Direct  computation  of  radiation  heat  generation  rates, 

2.  Greater  resolution  of  gamma-ray  energy  spectra. 

Capability  for  testing  more  complex  geometries. 
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The  logic  involved  in  the  programed  solution  of  the  spectral 
and  heat  and  dose  equations  as  well  as  the  equations  themselves 
are  described  in  Section  II,  and  the  instructions  and  data  formats 
required  in  order  to  use  these  programs  are  given  in  Section  III. 
Pour  appendices  contain* 

1.  A  list  of  the  symbols  used  in  the  flow  diagrams  for 
the  two  programs, 

2.  Derivations  for  the  neutron  flux-to-heat  conversion 
coefficients, 

3.  Tables  of  data  for  the  materials  libraries,  and 

4.  Neutron  r-eference-materlal  comparison. 
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II  PROaRAM  Loaic 


Each  penetration  program  Is  divided  Into  three  subprograms, 
namelys  geometry,  gamma,  and  neutron  routines.  The  gamma  And 
neutron  routines  are  the  same  In  both  programs.  The  geometry 
routine  of  the  first  code  (C17)  Is  restricted  to  geometries  com¬ 
posed  of  frustra  of  rectangular  pyramids  and  coaxial  cylinders  and 
their  annuli.  The  geometry  routine  of  the  second  program  (L-63) 
accepts  a  more  general  class  of  solids,  specifically,  cylinders  and 
their  annuli  which  are  defined  about  arbitrary  axes,  sectors  of 
these  cylinders,  and  frustra  of  pyramids  whose  bases  are  quadrilat¬ 
erals.  In  addition,  using  the  spherical  option,  spheres,  hollow 
spheres  or  hemispheres,  spherical  sectors,  and  spherical  sectors 
with  one  or  two  ends  cut  off  may  be  defined  for  L-63.  This  code 
also  accepts  regions  within  regions  and  regions  within  regions 
within  regions  In  which  the  geometry  types  can  be  varied.  These 
routines  are  described  In  the  following  sections. 

2.1  Geometry  Calculation 

The  Intercepted  distances  In  each  material  along  the  llne- 
of-slght  Joining  each  source  point  with  each  receiver  point  are 
required  In  the  gamma  and  neutron  routines.  The  purpose  of  the 
geometry  routine  Is  to  compute  this  Information  from  data  which 
describe  the  geometry  of  tne  system  and  the  location  of  source 
points  and  detector  pclnts. 

The  method  described  here  Is  not  the  classical  method  of 
determining  the  Intercepted  distances  between  two  points,  rather 
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it  is  an  iterative  method  based  upon  the  stepping-point  concept. 
Using  this  method,  a  point  is  moved  along  a  line  in  steps  of 
known  length.  The  stepping  point  P  (with  components  Xp,  y  ,  z  ) 
is  originally  coincident  with  the  source  point  S  (with  components 
Xg,  yg,  Zg).  A  stsp  slong  the  llne-of-slght  toward  the  detector 
point  is  accomplished  by  adding  a  constant  K  times  the  direction 
cosine  vector  L  to  the  stepping-point  vector.  Thus, 

T — y  T  -f  kIT, 

(The  above  expression  should  be  read  as;'*P  is  replaced  by  P'f'KL". ) 
Or  component  wise, 

Xp - *■  Xp  + 

yp — >  yp  + 

^p— 

where  ^3  direction  cosines  of  the  source  detector 

llne-of-slght  in  the  x,y,  and  z  directions,  respectively. 

The  stepping  point  must  be  identified  as  being  in  some 
particular  geometric  volume  before  the  first  step  is  taken.  Then, 
after  each  step,  the  stepping  point  must  be  tested  again  to  ascer¬ 
tain  whether  it  is  still  in  the  same  volume.  Once  a  boundary  has 
been  crossed,  a  vernier  effect  may  be  achieved  by  taking  one  step 
back,  reducing  the  step  size,  and  repeating  the  procedure  until 
the  difference  between  the  stepping -point  position  and  the  detector 
side  boundary  of  the  volume  is  less  than  the  boundary  uncertainty 
parameter, 
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The  method  of  defining  the  geometric  volumes  for  these  pro¬ 
grams  rests  on  the  concept  of  an  x-plane.  As  used  here,  for  a 
given  Cartesian  coordinate  system,  an  x-plane  Is  a  set  of  numbers 
sufficient  to  define  a  plane  area  perpendicular  to  the  x-axls  of 
the  system  plus  the  x-coordlnate  of  this  area,  or  the  areas, 
depending  upon  the  context.  Two  examples  of  x-planes  are  shown 
In  the  sketches  below. 
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The  first  of  these  shows  a  rectangular  area,  hence  the  x-plane 
consists  of  the  set  x,  Yraax*  Zmin'  ^max*  or,  the  rectangular 

area  at  x  bounded  by  the  last  four  of  the  above  numbers.  The 
second  Is  a  circular  annulus,  so  that  the  x-plane  consists  of  the 
set  X,  Rin»  or  the  annular  area  at  x  bounded  by  and  R^. 

For  these  program? ^  a  volume  Is  defined  by  two  or  more  x-planes 
of  the  same  type,  and  that  portion  of  the  volume  surface  which  Is 
not  coincident  with  one  of  the  defining  x-planes  Is  defined  from 
either  a  linear  Interpolation  or,  for  the  spherical  option, 
a  particular  second  order  interpolation  between  similar  points 
of  adjacent  x-planes.  A  volume  defined  by  rectangular  x-planes 
is  sketched  below.  This  volume  Is  defined  by  three  x-planes  of 
the  type  depicted  In  the  first  example  above. 
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The  class  of  "spherical  volumes"  acceptable  to  the  complex-geometry 
routine  consists  of  those  whose  centers  lie  on  the  x  axis  (see  sketch). 
Thus,  the  equation  for  this  class  Is 

(x-K)2  -h  =  0, 

where 

K  Is  the  x-coordlnate  of  the  center,  and 
Z’  Is  the  radius  of  the  sphere. 

The  Interpolation  formula  Is  derived  below. 


Assume  a  portion  of  a  sphere  Is  to  be  defined  with  center  at 
X(,,  radius  ,  and  truncated  at  x^^  and  x^.  Then  =^^^-(xcc 
and  -Xg  )^.  Then,  for  some  Xp,  such  that  x^f.  ^  x^fi  , 

the  point  R  may  be  found  In  terms  of  Xp,  3^  ,  R^  #  »  and  R^  by 

successive  application  of  the  Pythagorean  Theorem: 

r2  =/>2-(Xp  - 

4 
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Combining  the  above  equations  gives: 


p2 

A  a 


( Ro<  ^  ^  ^  ^  ^  ^  ^  ^oc  ~  ^ 

- 3J.___ - 


2 

X 


or 


A  volume  defined  by  cyllndrlc  x-planes  of  the  type  depicted  In  the 
second  example  is  shown  below.  It  should  be  noted  that  the  axis  of 
symmetry  of  this  figure  is  coincident  with  the  x  -axis  of  the 
coordinate  system  used  to  define  the  volume.  This  restriction  holds 
for  the  definition  of  all  volumes  used  In  the  programs,  excepting 
those  composed  of  x-planes  of  the  type  shown  In  the  first  example 
above  and  the  complex  Cartesian  x-planes  for  the  comp lex -geometry 
program.  Different  coordinate  systems  may  be  used  for  different 
sets  of  volumes  in  the  defined  region  by  defining  them  with  respect 
to  a  reference  system.  Z 


Examples  of  spherical  volumes  which  may  be  defined  by  a  coordinate 

type  (CT)  of  -2  or  -3  where  the  required  x-plane  parameters  arei 

1.  Sphere  with  center  at  x  =  a«  and  radius  b: 

CT  ■  -2  (simple  cylindric  x-planes) 

First  x-plane  t  x  r  a  -  b,  s  0,  R^  =  0 

Second  x-plane  t  x  =  a  ^  b,  Rj^n  =  0,  r”  =  0 

Y 


2.  Sector  of  a  hollow  sphere  or  spherical  shell  with  center  at 
a«  inner  radius  b^  and  outer  radius  t>2t  with  the  sector 
occupying  the  azimuthal  anglular  range  (in  the  yz-plane) 

^2,  to  02  (see  sketch  on  next  page): 

CT  s  -3  (complex  cylindric  x-planes) 

First  x-plane  t  x  s  a  -  b2*  R^^jj  r  0#  R©  =  ^  =  ^1'  ■  02 

Second  x-plane:  x  s  a  -  b^^,  Rj^p,  =  0,  R^  2  ^b|-b|,  =  0\&2~02 
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Third  x-plane* 
Fourth  x-planes 


X  ■  a+bj^,  0,  zvb2“b2^,  @2,  *  0\*  “  ^2 

X  s  a*b2#  R^n  z  0#  Rq  =  ~  ^1*  =  ^2 


Y' 


Note:  The  primed  coordinate  system  Is  shown  to  Illustrate  the 
method  of  measuring  the  azimuthal  angles  and  §2.  This  system 
Is  merely  a  translation  along  the  x  axis  of  the  unprlmed  system. 


For  clarity  and  conciseness,  the  following  terms  are  used  In  des¬ 
cribing  both  of  the  geometry  routines,  although  those  terms  pre¬ 
ceded  by  an  asterisk  are  used  for  the  complex  geometry  only.  The 
symbols  used  In  the  geometry  flow  diagrams  are  described  In  Appendix  A. 

1.  x-plane  -  this  concept  Is  defined  above. 

2.  element  -  a  volume  defined  by  two  adjacent  x-planes. 

3.  volume  -  a  set  of  adjacent  elements  defined  by  a  con¬ 
sistent  set  of  x-planes. 

Three  types  of  volumes  are  defined  for  the  complex-geometry 
routine  t 

a.  *subreglon  -  a  volume  containing  only  one  material 
and  containing  no  other  volume. 

b.  *reglon  -  a  volume  which  may  contain  an  arbitrary 
number  of  subregions  of  various  materials  and  has 
associated  with  It  a  base  material  which  occupies 
that  portion  of  the  region  not  occupied  by  sub- 
regions  • 

c.  ^master  region  -  a  volume  which  may  contain  up  to 
10  regions  of  various  kinds  and  has  associated  with 
It  a  base  material  similar  to  that  for  regions. 

Notet  All  regions  should  be  contained  In  a  master 

region,  and  all  subregions  should  be  contained 
In  a  region. 

For  the  simple-geometry  routine,  only  one  type  of  volume 
Is  defined,  and  It  Is  called  a  region.  For  this  routine  there 
Is  no  volume  containment. 

4.  dummy  volume  -  the  routine  requires  that  every  part 
of  the  problem  space  have  some  material  associated 
with  It.  The  dummy  volume  Is  not  a  defined  volume  In 
the  sense  that  a  number  of  x-planes  aj?e  required  to 
delineate  this  volume.  Rather,  It  Is  assumed  that 
the  dummy  volume  occupies  all  space  not  occupied  by 
other  defined  volumes.  For  the  complex-geometry 
routine,  the  last  master  region  defined  In  the 
geometric  Input  must  be  the  dummy  volume,  while 

for  the  simple-geometry  routine,  the  last  region 
Is  the  dummy  volume. 
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source~detector  line  -  the  straight  line  connecting  a 
source  point  and  a  detector.  (If  this  line  has  dlrec- 
tlon«  as  In  computing  the  direction  cosines  of  this 
llne«  the  sense  Is  from  source  to  detector.) 

6.  segment  -  the  penetration  distance  along  the  source- 
decector  line  through  one  volume  multiplied  by  the 
density  of  the  material  of  the  volume «  or  for  two  or 
more  volumes  of  the  same  material  which  are  either 
adjacent*  or  separated  only  by  void*  the  associated 
segment  Is  the  sum  of  the  products  of  the  penetration 
distance  through  each  volume  and  the  density  of  the 
material  therein.  (No  segment  or  portion  thereof 
corresponds  to  a  void  volume.) 


That  la*  If  ^  and  M.^n*  Mj,  arid  M 

are  non-zero  whers  Mt  Is  a  number  used  to  Identify 
a  material  In  the  itn  volume  along  a  source-detector 
line  and  «  0  If  the  material  Is  a  volc^  the  segment 
Wj^  Is  given  by 


Wi  -  ti/>i  * 

where  Is  the  density  of  the  material,  and  t  Is  the 
thickness  (along  the  source -detector  line). 

If  2*  1  “  ®  (void)*  then 

Wl  =  (HZ'  lv-2^- 

7.  detector-side  boundary*  source-side  bounda^  (of  a 
y^ume)  -  the  two  adjacent  Intersections  of  the  source- 
detector  line  and  the  volume  surface*  the  source-side 
boundary  being  closest  to  the  source  point. 

8.  *envelope  -  for  cartesian  volumes  -  the  smallest  rectan- 
gular  parallelepiped  which  contains  the  volume  In 
question.  This  envelope  Is  described  by  two  x-values* 
two  y-values*  and  two  z-values. 

for  cylindrical  volumes  -  the  smallest 
volume*  described  by  two  x-planes  of  the  same  type  used 
to  define  the  volume*  which  contains  the  volume. 

9.  *base  material  -  for  a  volume  -  that  material  which 
fills  the  volume*  except  for  the  space  occupied  by 
subvolumes. 


-  for  a  problem  -  that  material  In  which 
the  system  Is  assumed  to  be  Immersed .  (This  Is  the 
material  associated  with  the  dummy  volume.) 
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2.1.1  Slmple-Qeometry  Routine 


2. 1.1.1  Volumes  Acceptable  to  the  Routine.  Each  volume  used 
to  describe  the  system  being  studied  should  be  defined  by  a  set 

of  x-planes  of  one  and  only  one  of  the  following  types  (as  Illustrated 
In  Figure  2-1). 

1.  Simple  Cartesian.  The  defined  area  Is  a  rectangle 
whose  sides  are  parallel  to  the  y  and  z  axes.  The 
numbers  necessary  to  define  the  area  are  Xp, 

*««<  Vn' 

Restrictions t  Orientation  -  ^  ^max 

^In  ^  ^xoAx 

2.  Simple  Cylinder.  The  defined  area  Is  a  circular  annulus. 

The  numbers  necessary  to  define  the  area  are  Xp« 

and  R^  (Pig.  2-1). 

"  Restrictions!  Orientation  -  All  cylinders  must  be 

defined  to  be  coaxial  with  the  x-axls 

of  the  coordinate  system  used  to  define 

the  system  and  R^  ^  R  . 

in  o 

2. 1.1.2  Method  Outline.  The  segments  and  the  associated 
material  Identification  numbers  are  evaluated  using  the 
stepplng-polnt  described  In  Section  I.  The  overall  flow  dlaigram 
for  the  simple  geometry  code  (C-17)  Is  shown  In  Figure  2-2.  A 
flow  diagram  for  the  geometry  routine  Is  shown  In  Figure  2-3. 
Nomenclature  Is  given  In  Appendix  A. 

For  a  given  detector  point,  all  and  the  number  of 

segments.,  1^^^  are  determined  for  each  source  point  by  the  method 
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(increasing  x  is  out  of  the  paper) 


NPC  14,512 


Cartesian; 

a.  simple-  b*  complex- 


Cyllndrical; 

c.  simple  <i,  complex  e*  complex-rectilinear 


f.  Combination  I  g.  Combination  II 


FIOUM  2-1.  X-MANE  GROSS  SECTIONS 
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FIOURE  2-2.  SIMME- GEOMETRY  OVERALL  FLOW  MAGRAM 
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outlined  below.  At  this  point,  the  neutron  and/or  gamma  calcula¬ 
tions  are  made,  and  this  step  Is  repeated  until  the  calculations  are 
made  for  each  detector  point. 

1.  A  source  point  is  chosen;  then  the  source-detector 
distance,  D,  and  the  direction  cosines  of  the  line 
from  the  source  to  the  detector  are  computed,  and  the 
stepping  point  Is  taken  to  be  the  source  point. 

2.  The  stepping  point  Is  tested  against  each  region <  to 
determine  the  region  In  which  the  stepping  point  lies. 

3.  If  the  stepping  point  Is  found  to  lie  within  one  of 

the  defined  volumes,  the  detector-side  boundary  Is 
found  by  the  method  outlined  above  and  the  stepping 
point  Is  put  past  the  detector-side  boundary  of  this 
volume  but  within  of  the  boundary.  Then  the 

distance  between  the  source-side  and  detector-side 
boundaries  Is  taken  to  be  a  tentative  segment.  The 
material  number  of  this  volume  Is  found,  and  If  It 
Is  non-zero  (material  number  zero  Is  taken  to  be 
void),  the  material  number  Is  tested  against  the 
last  non-zero  material  number.  If  the  material 
numbers  are  equal,  the  materials  are  the  same,  and 
the  tentative  segment  Is  multiplied  by  the  appro¬ 
priate  density  and  added  to  the  last  segment.  If 
the  material  numbers  are  not  equal,  the  materials 
are  different,  and  the  tentative  segment  Is  multi¬ 
plied  by  the  appropriate  density  and  stored  as  a 
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segment  along  with  the  corresponding  material 
number  and  segment  counter  1.  If  the  material 
number  Is  zero,  the  tentative  segment  Is  not  used 
as  a  segment  or  a  portion  thereof. 

4.  Steps  2  and  3  are  repeated  until  either 

a.  the  stepping  point  Is  found  to  have  passed  the 
detector  -  In  this  case,  the  distance  between 
the  source-side  boundary  and  the  detector  Is 
called  a  tentative  segment,  which  Is  treated 
as  In  Step  3  above,  and  the  program  goes  to 
Step  1,  unless  the  source  points  have  been 
exhausted,  or 

b.  the  stepping  point  Is  found  to  be  In  an  unde¬ 
fined  region  (the  dummy  material)  by  testing 
against  all  defined  regions  (Step  2)  and  find¬ 
ing  the  stepping  point  In  none  of  them  -  In 
this  case,  the  stepping  point  Is  taken  to  be 
at  the  detector,  and  the  distance  between  the 
last  detector-side  boundary,  and  the  dectector 
Is  taken  ^'o  be  a  tentative  segment  and  Is 
treated  as  In  Step  3>  Then  the  program  goes 
back  to  Step  1  unless  the  source  points  have 
been  exhausted. 

2.1.2  Complex-Geometry  Routine 

2. 1.2.1  Geometric  Ordering  of  Volumes.  In  order  to  facili¬ 
tate  the  region-search  routine,  each  volume  should  be  contained 
in  one  of  ten  master  regions.  Further,  If  geometric  complexity 
warrants,  each  master  region  may  contain  up  to  ten  regions  and 
each  region  an  arbitrary  number  of  subregions.  This  mode  of 
categorization  has  the  advantage  of  reducing  the  number  of 
regions  the  routine  must  search  at  any  particular  boundary. 

Also,  each  master  region  and  each  region  will  have  associated 
with  It  a  base  material  (which  should  be  either  the  predominant 
material  therein  or  that  material  which  would  minimize  the 
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geometry  description),  so  that  only  those  volumes  containing 
materials  other  than  the  base  need  be  defined. 


2. 1.2.2  Volumes  Acceptable  to  the  Routine.  Each  volume 
used  to  describe  the  system  being  studied  should  be  defined  by 
a  set  of  x-planes  of  one  and  only  one  of  the  following  types 
Illustrated  In  Figure  2-1. 

1.  Cartesian 

a.  Simple.  The  defined  area  Is  a  rectangle  whose 
sides  are  parallel  to  the  y  and  z  axis.  The 
numbers  necessary  to  define  the  area  aret  x  , 

^mln'  ^max*  ^mln*  ^max  2-la).  ^ 

Restrictions:  Orientation  -  and 

niin  —  tnax 

^mln^^max* 

b.  Complex.  The  defined  area  Is  a  quadrilateral. 

The  numbers  necessary  to  define  the  area  are: 

Xp,  Y^,  Z^,  Yg,  Zg,  Yy  Zy  Y^,  Z^^  (Pig.  2-lb). 

Restrictions:  On  the  points  s  (Xp,Y^,Z^), 

1  =  1, 2,3,4 

Orientation:  P^^Pg  and  P^P^^  are  not 
parallel  to  the  z  axis. 

^1^4  and  parallel  to 

the  y  axis. 

mln(Z2^Zg)^  mln(Z^Z2^),  and 

mln(Y^Y2|)^  mlnCYgY^)* 

Connectivity:  P,  Is  adjacent  to  P2 
and  P^^  ^ 

Convexity:  All  Internal  angles  less 
than  1800. 
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2.  Cylindrical 


a.  Simple.  The  defined  area  Is  a  circular  annulus. 

The  numbers  necessary  to  define  the  area  are  x  , 

and  (Pig.  2-lc).  P 

b.  Complex.  The  defined  area  Is  a  sector  of  a  cir¬ 

cular  annulus.  The  numbers  necessary  to  define 
the  area  are  Xp,  and  R^  (Plg.  2-ld). 

c.  Complex  rectilinear.  The  defined  area  Is  a 
quadrilateral,  two  of  whose  sides  are  on  radii 

of  a  circle.  The  numbers  necessary  to  define  the 
area  ares 

Xp,  @2,  Rj  .  Ri  ,  Rp^,  and  Rq^  (Pig*  2-le). 

d.  Combination.  The  defined  area  Is  a  generalized 
quadrilateral,  two  of  whose  sides  are  on  radii  of 
a  circle.  The  third  side  Is  a  circular  arc,  and 
the  fourth,  a  straight  line. 

I.  The  arc  forms  the  Inner  side.  The  numbers 
necessary  to  define  the  area  are :  x  , 

Rln»  Roi*  and  R^^  (Pig*  2-lf).  V  i  ^ 

II.  The  arc  forms  the  outer  side.  The  numbers 
necessary  to  define  the  area  arei 

**lni'  %n2»  **o  2-lg). 

The  defined  volumes  using  these  x-planes  may  be  defined  In  terms 

of  any  arbltraz^^  coordinate  system.  The  only  restriction  Is  that 

the  7  axis  of  a  coordinate  system  used  to  define  a  cylindrical 

volume  must  be  the  axis  of  the  cylinder. 

Restrictions:  Orientation  -  All  cylinders  must 
be  defined  to  be  coaxial  with  the 
X  axis  of  the  coordinate  system 
used  to  define  the  cylinder.  All 
azimuthal  angles  ar^  positive  and 
measured  from  the  y  axis. 

Rlni-  "oi*  ^ 


All  angles  are  to  given  In  degrees. 
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2, 1,2. 3  Method  Outline.  The  segments  and  the  associated 
material  Identification  numbers  are  evaluated  using  a  modifica¬ 
tion  of  the  stepplng-polnt  method  described  above.  Plow  diagrams 
of  the  geometry  routine  are  shown  In  Figures  2-4  through  2-9. 

The  general  method  of  evaluating  Wj^  and  Is  as  follows: 

1.  The  envelopes  for  all  volumes  are  computed  and  the 
coordinate  transformation  matrices  (to  transform  a 
point  from  the  base  system  to  the  given  system)  are 
computed  for  all  coordinate  systems.  (Note:  since 
the  matrix  A  represents  a  rotational  trans format Ion « 

It  Is  orthogonal,  l.e..  Its  Inverse,  equals  Its 

transpose.  A*’’.)  This  procedure  Is  shown  In  Figure  2-4. 

2.  A  detector  point  Is  chosen  and  Its  coordinates  are 
transformed  Into  Cartesian  coordinates  and  Into  the 
base  system  as  shown  In  Figure  2-3,  and  all  Wj^,  Mj^, 
and  lina-r  determined  for  every  source  point  by 
the  method  outlined  In  Steps  3f  4,  and  3*  At  this 
point,  the  neutron  and/or  gamma  calculations  are 
made  and  this  step  Is  repeated  until  the  calculations 
are  made  for  each  dete^'.tor  point. 

3.  A  source  point  Is  chosen  and  Its  coordinates  are  trans¬ 
formed  Into  Cartesian  coordinates  and  Into  the  base 
system.  Then,  the  source-detector  distance  and 
direction  cosines  are  computed.  The  stepping  point 

Is  taken  to  be  the  source  point  (Fig.  2-6). 
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4,  The  stepping  point  is  tested  against  each  master  region 
to  find  the  master  region  in  which  the  stepping  point  lies 
(Pigs,  2-7  and  2-9).  The  stepping  point  is  transformed 
into  the  coordinate  system  of  the  volume  and  tested 
against  the  envelope.  If  the  point  is  not  within  the 
envelope^  a  new  volume  is  chosen  and  the  procedure  begins 
over.  If  the  point  is  within  the  envelope^  then  the  point 
is  tested  to  see  if  it  is  in  the  volume  (Fig.  2-9).  If 
the  stepping  point  is  not  in  the  volume,  a  new  volume 
is  chosen  and  tested.  If  the  stepping  point  is  in  the 
volume,  the  detector-side  boundary  of  this  volume  is 
found  by  the  stepping-point  method,  and  the  distance  t 
between  this  point  and  the  last  boundary  point  is 
determined  (for  the  first  boundary,  the  source  point  is 
used).  The  volume  is  then  tested  to  see  if  it  contains 
subvolumes.  If  it  does  not,  t  is  multiplied  by  the 
material  density  associated  with  the  volume  and  stored 
as  the  1^^  segment  W^;  Mj^,  the  material  identification 
number,  is  also  stored.  Then,  with  the  stepping  point 
barely  beyond  the  last  computed  boundary  point,  this 
step  is  repeated.  If  the  volume  is  found  to  contain  a 
subvolume  (Fig.  2-8),  the  stepping  point  is  put  on  the 
source-side  boundary  point  of  the  volume,  and  if  the 
containing  volume  is  a  master  region,  this  step  is 
repeated,  testing  against  the  regions  contained  within 
this  master  region;  or  if  the  containing  volume  is  a 

region,  this  step  is  repeated,  testing  against  the 

36 


subregions  contained  In  this  region. 

3.  After  the  last  step  has  been  performed  a  number  of 
times,  one  of  two  things  will  occur:  Upon  testing 
the  volumes  to  find  which  volume  the  stepping  point 
Is  In,  It  will  be  found,  by  elimination,  that  either 
the  stepping  point  Is  In  none  of  the  volumes,  or.  If 
the  set  of  volumes  are  master  regions.  It  may  be 
found  that  the  stepping  point  has  gone  past  the 
detector  point. 

In  the  first  Instance  (Pig.  2-8),  the  stepping  point 
Is  assumed  to  be  In  the  base  material  and  the  portion 
of  the  source-detector  line  lying  between  the  stepping 
point  and  the  detector-side  boundary  point  of  the  con¬ 
taining  region  (or  the  detector  point  If  the  stepping 
point  Is  outside  all  master  regions)  Is  tested  to  see 
If  this  line  passes  through  a  volume  (or  volumes). 

If  It  does  not,  the  stepping  point  Is  advanced  to  the 
detector-side  boundary  of  the  containing  regions  and 
the  distance  through  the  base  suiterlal  Is  taken  to  be 
a  segment  Wj^,  and  and  the  associated  Mj^  ax«  stored. 
If  the  volumes  of  the  set  being  tested  are  the  master 
]?eglons,  Igj  Is  stored,  and  the  program  goes  to  Step  3* 
If  the  volumes  of  the  set  being  tested  are  regions 
(subregions),  then  with  the  stepping  point  Just  past 
the  detector-side  boundary  of  the  containing  volume. 
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the  set  of  master  regions  (regions)  Is  considered <  and 
the  program  goes  to  Step  4, 

If  a  volume  (or  volumes)  lie  on  this  portion  of  the 
source-detector  line,  then  the  source-side  boundary 
point  of  the  volume  closest  to  the  source  Is  found, 
and  the  distance  that  the  stepping  point  goes  through 
the  base  material  Is  called  a  segment  Vl±,  Wj^  and  the 
associated  are  stored.  The  program  then  goes  back 
to  Step  4. 

If  the  stepping  point  has  gone  past  the  detector,  the 
segment  Is  taken  to  be  the  portion  between  the 
volume  source-side  boundary  point  and  the  detector. 

The  program  then  goes  back  to  Step  3« 


38 


q  -*  0, 
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FIGURE  2-4.  ENVELOPE  PARAMETERS  AND  COOR 
TRANSPORTATION  ROUTINE 


Computation  of  the  U*'3  and  V*-'s, 
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(URE  2*4.  ENVELOPE  PARAMETERS  AND  COORDINATE 
TRANSPORTATION  ROUTINE 


Find  Envelope 
Parameters  and 
Transformation 
Matrices 
(Fig.  2-4) 


IPLEX- GEOMETRY  OVERALL  FLOW  DIAGRAM 
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FIGURE  2-7.  ROUNOARY  DETERMINATION  SURROUTINI 


BOUNDARY  DETERMINATION  SUBROUTINE 
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(to  Pig. 2-7) 


FIGURE  2-9.  ELEMENT  DETERMINATION 


2.2  Calculation  of  the  Qymna-Ray  Number-Flux  Bnergy  Spectrum. 

TSbse  ftate^  and  ieat  fleneration  ftate 

The  gamma-ray  number-flux  energy  spectrum  at  a  point  In  space 
may  be  divided  Into  two  components,  a  direct-beam  portion  and  a 
scattered  portion.  Thus,  F(Eg^),  the  spectral  point  for  energy 
at  a  point  Isotropic  detector  located  at  a  point  r  In  an 
Infinite  medium  due  to  photons  emitted  from  a  point  Isotropic 
source  located  at  the  origin,  and  with  source  spectrum  S(E|,), 

Is  given  by 

P(Ea)  *  -  +1  S(Eb)  l'(Eb,Ea,r)  dE^,  (l) 

where  the  first  term  on  the  right-hand  side  of  the  equation  repre¬ 
sents  the  direct-beam  component,  and/4.(Ea)  Is  the  linear  attenua¬ 
tion  coefficient  for  photons  with  energy  E^  In  the  material  In 
question.  The  second  term  In  Equation  1  represents  the  scattered 
portion  of  the  spectrum;  the  function  1  ^(Eb«Ea»r)  Is  the 
scattered  portion  of  the  differential  number  spectra  of  the  photon 
number  flux  In  an  lnfl.nlte  medium.  The  Integration  In  the  second 
term  Is  over  all  Initial  photon  energies  greater  than  or  equal  to 
Ea» 

The  differential  energy  spectra  of  the  gamma-ray  number 
flux  at  a  point  Isotropic  detector  resulting  from  a  point  Iso¬ 
tropic  source I  In  an  Infinite  medium,  as  computed  by  the  moments 
method  (Ref.  2)^  have  been  used  In  determining  the  function 

l'(Eb,Ea#r).  These  differential  energy  spectra  Iq  may  be  defined 
as  probabilities,  per  unit  Initial  photon  energy,  that  a  quantum 
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of  the  gamma-ray  energy  flux  will  be  degraded  from  an  Initial 


energy  to  an  energy  time  the  quantum 

reaches  a  point  isotropic  detector  at  r.  Hence,  the  differential 
energy  spectra  and  the  differential  number  spectra  are  related 
by 

- i" - =  I  (Bb-V'’)- 


The  data  tabulated  In  Reference  2  are  in  the  form 

f  z  4frr^  e^(Eb)r  i^(Ej^,Eg^,i((E^j)r)  (3) 

for  discrete  values  of  the  three  variables  and  for  various  materials. 
Thus,  when  the  Integral  Is  approximated  by  a  finite  sum,  the 
equation  to  be  solved,  from  Equations  1  and  2  and  the  definition 
above  (Eq.  3)»  Is 


-A<(E  )r  _  -A(Eb)r 


(4) 


krt  r‘ 


4Tt'r*^  E, 


a 


where  hf(E_,E.  )  is  a  numerical  integrating,  or  histogram,  factor, 
a  D 

and  the  summation  Is  over  Initial  photon  energy  E^^,  and  is  from 
the  energy  E^^  to  the  maximum  value  of  Initial  photon  energy. 

Equation  4  defines  the  energy  spectrum  of  the  gamma-ray 
number  flux  In  an  Infinite  medium  In  terms  of  the  spectrum  at  a 
point  detector  at  a  distance  r  from  the  point  Isotropic  source. 
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In  order  to  apply  this  method  to  determining  spectra  In  and  around 
reactor  shield  systems.  It  Is  necessary  to  assume  that  the  spectra 
may  be  reconstituted  at  each  boundary,  and  that  the  spectrum  on 
the  detector  side  of  a  boundary  Is  equal  to  the  spectrum  on  the 
source  side  of  the  boundary.^  As  an  example  of  this  Idea,  con¬ 
sider  a  source  and  detector  shown  In  the  sketch. 


Material  1 


Material  2 


Photons  from  the  source  S,  whose  spectrum  Is  taken  to  be  S2(E5), 
penetrate  the  two  slabs  with  thicknesses  tj^  and  t2  and  are  detected 
at  the  receiver  at  D.  The  spectrum  at  D,  as  an  application  of 
Equation  4,1s  given  byi 


e-^g(gb)t2 


Where  D  is  the  source -detector  distance,  and 

82  la  the  spectrum  (excluding  geometric  attenuation) 
at  the  boundary  between  the  two  materials. 

82  Is  given  In  terms  of  8^,  the  source  spectrum,  by 


S2(Ba)  =  Si(Ea)  e  3;^(Eb)  f 

Ec 


hf(Ea,E5) 


The  validity  of  this  assumption,  and  the  second  assumption 
that  follows,  has  been  tested  experimentally;  the  experi¬ 
ment  and  analysis  are  reported  In  References  5  and  6 . 
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A  second  assumption  inherent  in  this  method  is  that  differential 
energy  spectra  for  infinite  media  may  be  used  to  describe  radiation 
transport  in  finite  media.  One  way  of  estimating  the  Importance  of 
this  flnlteness  of  geometry^  is  to  assume  the  correction  to  be 
independent  of  final  energy  and  to  multiply  the  differential  energy 
spectra  by  the  edge  correction  factors  of  Berger  and  Doggett 
(Ref.  7). 

These  factors  are  defined  as 


g 


r)/Ea^dEa 
r)/EgJ  dEa 


(5) 


where  Io(E^j,Ej^,r)  and  i^(E|j,Ej^,r)  are  differential  energy  spectra 
for  a  finite  slab  and  for  an  infinite  medium,  respectively.  The 
function  g  may  be  rewritten  as 


0 


It  is  assumed  that  the  ratio  Iq/Iq  is  Independent  of  degraded 
energy  Ej^  and,  hence,  that  the  edge  corrections  are  Just  the 
ratio  of  the  differential  energy  spectra  for  a  finite  medium 
and  the  spectra  for  an  infinite  medium: 


Thus,  if  edge  corrections  are  to  be  used,  the  function  f  In  Equation 
4  should  be  replaced  by  fg.  This  correction  is  optional  in  the 
codes . 

It  has  been  assumed  that  any  source  may  be  replaced  by  a 
set  of  point  isotropic  sources^  hence,  the  spectrum  from  the  total 
source  is  taken  to  be  the  sum,  over  all  source  points,  of  the 
spectra  due  to  each  point.  The  points  of  this  total  spectrum  are 
then  multiplied  by  either  the  flux-to-dose  or  flux-to-heat  conver¬ 
sion  factors  and  by  the  appropriate  integrating  or  histogram  factors, 
and  summed,  over  energy,  to  give  either  dose  rate  or  heat  genera¬ 
tion  rate. 

2.2.1  Spectral  Calculation  (for  one  source  point) 

The  method  used  to  evaluate  the  spectral  equation  (£q.  4)  is 
shown  in  Figure  2-10.  This  flow  diagram  shows  the  method  for 
obtaining  the  spectrum  at  a  detector  from  one  source  point,  using 
the  segments  determined  by  the  geometry  routine. 

The  evaluation  of  the  spectral  points  ^  is  as  follows i 

1.  The  spectrum  associated  with  the  source  polnt^  (from 
Library  3)  is  taken  to  be  the  working  spectrum  St,  and 
the  first  segment  and  the  material  number  for  this 
segment  are  set  up  to  be  used.  Note:  the  subscript  1, 
on  W  and  M  is  used  to  denote  segment  number  . 

2.  The  subscript  a,  which  denotes  final  energy  (and,  in  a 
manner  similar  to  b,  increases  as  final  energy  decreases), 

1  The  subscript  b  denotes  initial  energy,  and  b  Increases  as 
Initial  energy  decreases;  l.e.,  b  =  1  denotes  10  Mev,  b  s  2. 
denotes  9  Mev,  etc. 
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Is  set  equal  to  one.  the  scattered  term  of  the 

spectrum  Is  set  equal  to  zero  In  order  to  Initialize 
the  summation  over  Initial  energy  (Eq.  4). 

3.  The  value  of  a  Is  tested.  If  a  Is  less  than  13«  b  Is 
set  equal  to  a.  If  a  Is  greater  than  or  equal  to  13« 
b  Is  set  equal  to  13. 

4,  The  values  for  the  mass  attenuation  coefficient,/^,  and 

effective  atomic  number,  are  taken  from  Library  1 

for  b  and  M^.  comes  from  the  geometry  routine, and 
the  number  of  relaxation  lengths  of  material  penetrated, 
y'bi'i.  calculated.  The  coefficients  for  the 
differential  energy  sncctral  function^ 

are  found  (from  Library  2)  as  a  function  of  a  and  b,  and 
*/<b^l 

the  function  fe/  Is  evaluated  using  the  variables 
^eff  *“^/^ 

3.  The  parameter  B  Is  tested  to  see  If  edge  corrections 
are  required  for  this  calculation.  If  B  Is  greater 

2 

than  zero,  the  coefficients  for  the  edge  corrections 


^  The  differential  energy  spectra  have  been  fitted  to  a  quadratic 
InAijV  and  Z^ff  for  each  possible  combination  of  a  and  b  (since 
degraded  energy,  for  a  particular  case,  must  be  less  than  or,  at 
most,  equal  to  the  Initial  energy,  a^b).  This  representation 
has  the  fora  _  t  j 

It  has  been  found  necessary  to  section  the  curve  fits  In  Z.£.f 
Into  three  ranges (0-26,  26-74,  and  74-92)and  to  use  different 
sets  of  coefficients  for  each  range.  See  Reference  3* 

^  The  edge  corrections  have  been  fitted  to  a  quadratic  of  the  form 

for  each  value  of  b.  It  has  been  found  necessary  to  section  the 
curve  fits  In/^bV  Into  two  ranges  (0-4  and  4-20)  and  to  use 
different  coefficients  for  each  range. 
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6. 


are  found  (from  Library  2)  as  a  function  of  b,  and  the 
function  g  Is  evaluated  using  the  coefficients  Zgff 
and  y^bWj^.  Then, 


A  > 


Is  formed, 
given  by 

J 


If  B  Is  less  than  or  equal  to  zero. 


J  Is 


The  histogram  factor  and  final  energy  Eg^  are  found 

(from  Library  2)  and 


(6) 

cl 

Is  formed,  where  Is  the  current  value  of  the 

spectral  point  for  b.  Expression  6  has  the  form  of 
the  summation  of  Equation  4,  except  that  the  limits  on 
the  summation  are  from  Eg^  to  the  current  value  of  E^* 

7.  The  parameter  b  Is  tested,  and  If  b  is  greater  than  one, 
the  current  value  of  b  Is  decremented  by  one  and  the 
routine  goes  to  Step  4.  If  the  current  value  of  b  Is 
one,  Is  now  the  scattered  term  of  Equation  4  and 
the  direct-beam  term  Is  to  be  computed  next.  First, 
the  mass  attenuation  coefflclenty^ls  found  (from 
Library  1)  for  a  and  Then<^i«r^  Is  computed  and  the 

new  current  value  of  the  spectral  point  for  a,  S*  . 

a.,  1, 

Is  given  by 


■h 


1 


e 


6o 


Next,  a  Is  tested  and  if  a  Is  less  than  14«  a  Is  Incre¬ 
mented  by  one,  and  the  routine  goes  to  Step  3.  If  a  Is 
equal  to  14,  1  is  tested,  and  if  1  is  less  than 
(which  is  computed  by  the  geometry  routine),  1  is  incre¬ 
mented  by  one,  and  the  routine  goes  to  Step  2;  if  i  is 
equal  to  1^^^  the  function  ^is  computed  from 


where  I  is  the  source  Intensity  (from  Library  3)* 

D  is  the  source-detector  distance  (from  the 
geometry  routine),  and 

S*,lmax  iB  the  new  current  value  of  the  spectral  point 

for  a  and  for  the  last,  or  1^^^^^,  segment. 

The  function  p|  Is  computed  for  each  value  of  a,  and 

Is  the  evaluation  of  the  photon  number-flux  spectrum 

at  the  detector  due  to  the  given  source  point. 

2.2.2  Determination  of  the  Spectrum,  Dose  Rate,  and  Heat 
generation  Rate 

The  method  used  to  determine  .the  spectrum,  dose  rate,  and 
heat  generation  rate  Is  shown  In  Figure  2-11.  This  flow  diagram 
shows  the  method  for  obtaining  these  parameters  from  the  spectra 
determined  by  the  spectral  routine  (Sec.  2.2.1). 

The  method  Is  as  follows: 

1.  The  following  parameters  are  set  to  zero  In  order  to 

Initialize  the  various  summations:  b^*,  H*^*,  and  — 

'  a 

the  dose  rate,  heat  generation  rate,  and  spectrum. 
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respectively —summed  over  all  source  groups.  The  source 
group  index  S*  is  set  to  one. 

2.  The  spectrum  summed  over  all  source  points  in  source 
group  S*,  is  set  to  zero  in  order  to  Initialize  the 
summation  over  source  points.  The  source-point  index  S 
is  set  to  one. 

y 

3.  The  spectral  function  P®  for  source-point  S  of  source- 

cl 

group  S*  is  found  (Sec.  2.2.1). 

4.  The  degraded-energy  index  a  is  set  to  one;  then  Q,  the 
source-point  print  option  is  tested.  If  Q  is  one,  dose 
rates  and/or  heating  rates  due  to  each  source  point  will 
be  calculated  to  be  printed  out  in  addition  to  the  normal 
printout.  If  Q  is  two,  only  the  normal  printout  will  be 
made,  and  for  each  value  of  a,  P^  obtained  in  Step  3 
will  be  added  to  the  current  value  of  a][.  If  Q  is  one, 
the  dose  rate  due  to  each  source  point  will  be  calculated 
if  X~2  and  the  heat  generation  rate  will  be  calculated 

V 

if  X— 2.  In  addition,  for  each  value  of  a,  P^, obtained 
in  Step  3«  will  be  added  to  the  current  value  of  a][. 

The  equations  used  to  compute  the  dose  rates  and  heat 
generation  rates  aret 

.nd  //y  . 

where  h^  are  integrating  or  histogram  factors  from 

Library  2, 

are  flux-to-dose  conversion  factors  from 
Library  2, 
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*re  flux-to-heat  conversion  factors  for 

material  of  unit  density,  and  are  a  function 

of  (from  the  geometry  routine)  as  well 

imax 

as  a^  they  are  listed  In  Library  1,  and 
Is  the  material  density  associated  with 
from  the  problem  deck. 


max 


3.  S,  the  source  point  Index,  Is  tested;  If  S  Is  less  than 

Sms#  total  number  of  source  points  In  the  source 

group,  S  Is  Incremented  by  one  and  the  routine  goes  to 

Y 

Step  3.  If  S  Is  greater  than  or  equal  to  S„„*  0-  Is 

ins  ^  cl 

now  the  total  spectrum  for  group  S*  and  the  dose  rates 
and/or  heat  generation  rate  are  calculated  as  In  Step  4 
with  0^  substituted  for  P^.  If  X  — 2,  D^*  Is  replaced  by 
the  sum  of  the  current  value  of  D^*  and  the  dose  rate 

due  to  the  s*^^  source  group.  If  X^2,  Is  replaced 
by  the  sum  of  the  current  value  of  and  the  heat 
generation  rates  due  to  the  S*^^  source  group.  In  any 
case,  0^  Is  replaced  by  the  sum  of  the  current  value 
of  and  0^  In  order  to  generate  the  total  spectrum. 

S*  Is  then  tested,  and  If  S*  Is  less  than  S*m«  the  total 
number  of  source  groups,  S*  Is  Incremented  by  one  and 
the  routine  goes  to  Step  2.  If  S*  Is  greater  than  or 
equal  to  S^,  then  the  current  values  of  0**.  d''  *,  and 
are  the  total  spectrum,  dose  rate,  and  heat 
generation  rate,  respectively,  at  a  detector  due  to 
all  source  points. 
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The  data  generated  by  this  routine  for  pz*ograin  output  aret 

a.  0^  and  0^  (for  a  =  1  through  14)  -  the  spectra 
of  the  gamma-ray  number  flux  summed  over  each 
source  group  and  over  all  source  groups «  respect¬ 
ively. 

b.  and  D^*  -  these  are  generated  only  If  X  ^  2; 
they  are«  respectively,  the  dose  rate  due  to  each 
source  group,  and  the  total  dose  rate.  If  Q  Is 
one.  In  addition  to  the  above,  the  dose  rate  due 
to  each  source  point  Is  printed  out. 

c.  and  -  these  are  generated  only  If  X  >  2: 
they  are,  respectively,  the  heat  generation  due 
to  each  source  group,  and  the  total  heat  genera¬ 
tion  rate.  If  Q  Is  one.  In  addition  to  the  above, 
the  heat  generation  rate  due  to  each  source  point 
Is  printed  out. 
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2.3  Caloulfttlon  of  the  Neutron  Number-Flux  Energy  Spectrum.  Dose 
fette.  and  rieat  Generation  Rate' 

The  neutron  calculation  la  baaed  on  the  aaaumptlon  that  the 
momenta-method  doae  ratea  and  energy  apectra  of  the  faat-neutron 
number  flux  due  to  a  point  Isotropic  Watt  flaalon  aource  In  an 
Infinite  medium  of  a  few  aelected  reference  materlala  may  be  uaed 
to  determine  the  apectra  and  dose  rates  In  arbitrary  combinations 
of  arbitrary  materials. 

The  method  used  to  determine  these  fluxes  and  dose  ratea  la 
based  upon  the  concept  of  an  equivalent  segment  The  equiva¬ 

lent  segment  for  a  material  whose  actual  segment  Is  W  and  whose 
removal  cross  section  (per  unit  density)  Is  Is  given  by 


5:  R 

-= - W  , 

^  RR 


(1) 


where  ^RR  the  removal  cross  section  (per  unit  density)  for  the 
chosen  reference  material.  The  total  equivalent  segment  Is  the 
sum  of  the  equivalent  segments  for  each  segment  along  the  source- 
detector  line.  This  total  equivalent  segment  is  used  to  deter¬ 
mine  the  attenuation  factors  for  the  spectral  points  and  dose 
rates,  and  these  factors,  when  multiplied  by  the  source  Intensity, 
and  the  geometric  attenuation  factor  are  taken  to  be  the  values 
of  the  spectral  points  and  dose  rate  at  the  chosen  detector  due 
to  one  source  point.  Hence,  the  spectral  points  and  dose  rate 
at  the  detector  due  to  the  total  source  is  the  sum,  over  all 
source  pointer  of  the  spectral-point  or  dose-rate  function  per 
source  point.  The  heat  generation  rate  is  computed  from  the  total 
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spectral  function  by  multiplying  It  by  appropriate  Integrating  or 
histogram  factors  and  the  fast-neutron  flux-to-heat  conversion 
factors  (Appendix  B)  and  summing  over  neutron  energy.  A  comparison 
of  the  various  reference  materials  In  a  given  geometry  Is  located 
in  Appendix  C. 

The  method  used  to  perform  the  neutron  calculations  Is  shown 
In  Figure  2-12.  This  flow  diagram  shows  the  method  for  obtaining 
these  parameters  from  moments-method  data  for  neutrons. 

The  method  used  Is  as  follows: 

1.  The  following  parameters  are  set  to  zero  In  order  to 

initialize  the  various  summations:  D^*  (for  a  =  1  to 
a^ax+-l)  ■  first  of  these  are  the  neutron 

spectral  terms,  and  the  term  is  the  dose-rate 

function  summed  over  source  points;  a  Is  the  neutron- 
energy  index  and  increases  with  decreasing  energy, 
i.e.,  a  s  1  for  neutron  heat  generation 

rate  summed  over  all  source  points.  S*,  the  source- 
group  Index,  Is  set  to  one. 

2.  The  following  parameters  are  set  to  zero  In  order  to 

Initialize  the  various  summations:  (for  a  =  1  to 

^max+1^  the  neutron  spectral  terms  and  dose  rate  and 
h"^,  the  heat  generation  rate  summed  over  all  source 
points  in  source  groups  S*.  S,  the  source  point  Index 
Is  set  to  one. 

3.  The  total  equivalent  segment  for  source  point  S  of 
source  group  S*  Is  found  from 


W 


eq 


■  ^  1  .  ^max  T  W. 

^‘4=11  ' 
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where  Is  the  number  of  segments  on  the  source- 

detector  line  (from  the  geometry  routine)^ 
\l±  Is  the  1^^  segment  as  computed  by  the 
geometry  routine^ 

^  Is  the  removal  cross-section  (per  unit 

«!■  V  «-h 

density)  for  the  1''"  material  and  Is  from 

Library  1  for  material  and 

^  Is  the  removal  cross  section  (per  unit 
*"RR 

density)  for  the  reference  material  and 
Is  from  Library  1  for  the  material  corre¬ 
sponding  to  the  reference  material. 

The  energy  Index  a  Is  set  to  one,  and  the  heat 

generation  rate  per  source  point.  Is  set  to  zero. 

4.  The  neutron  moments -met  hod  energy  spectra^  Pj^*  Is 
generated  and  the  spectral  point  for  energy  a. 


^  =  2.46  •  •  I  , 

Is  computed  (1  Is  the  source-point  fission  Intensity  from 
Library  3  and  2.46  Is  the  number  of  neutrons  per  fission). 
Then,  the  function  Is  replaced  by  the  current  value 

of  plus  P^  and  the  detector  print  option  X  Is  tested. 


^  The  moments -met  hod  attenuations  for  neut^ns  have  been  fitted  to 
a  polynomial  of  the  form  ^  =  expfCiW^-^  C2W3*f  CoW^-f 

Ci|W-«‘C3^and  the  coefficients  for  this  polynomial  are  listed 
In  Library  1  In  order  of  descending  energy;  the  last  set  Is 
for  the  dose-rate  curve.  It  has  been  found  necessary  to  section 
these  curve  fits  In  order  to  obtain  a  better  fltj  hence,  two 
sets  of  coefficients  are  shown  In  Library  1,  one  for  W*BP,  and 
one  for  w  2BP  (BP  Is  an  arbitrarily  chosen  segment,  and  for  the 
data  shown  In  Appendix  C  the  most  common  value  Is  60  gm/crn^. 
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If  X  Is  greater  than  or  equal  to  2,  the  source-point 
option  Q  Is  tested.  If  Q  Is  1,  then  h*'**  Is  replaced 
by  the  current  value  of  H"  plus  the  product  of  Pg, 
§”^(the  neutron  flux-to-heat  conversion  factors  from 
Library  1)«  and  hf^^  the  Integrating  or  histogram  factors 
from  Library  1. 

3.  a  Is  tested;  If  a  Is  less  than  routine  goes  to 

Step  4,  If  a  =  *niax*  ^  1>  and  X2  2,  the  current  value 

of  h”**  Is  the  heat  generation  rate  due  to  the  one  source 
point  and  this  value  Is  to  be  printed  out  by  the  program. 
If  a  =  ^niax  2,  the  dose  rate  Is  to  be  computed; 

hence,  a  Is  set  equal  to  &nd  Step  4  Is  repeated; 

then  Q  Is  tested,  and  If  Q  Is  1  the  dose  rate  for  the 
source  point,  printed  out  by  the  program. 

6.  The  source-point  Index  S  Is  tested.  If  S  is  less  than 
Smax(S*),  the  number  of  source  points  In  the  S*^*’  source 
group  (from  Library  3)«  then  S  Is  Incremented  by  one  and 
the  routine  goes  to  Step  3*  If  S  =  S^^y  (S*),  the 
current  values  of  (for  a  =  1  through  a^^^)  are  the 
spectral  points  for  the  energy  spectra  of  the  neutron 
number  flux  due  to  the  S*^^  source  group  and  the  current 
value  of  Is  the  dose  rate  for  this  group. 

These  numbers  are  printed  out  by  the  program.  If  X  ~  2, 
and  Q>1  the  heat  generation  rate  due  to  this  source 
group  Is  calculated  from 
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a 


max 


n 


=  T.  >^1  HJV. 


a -1 


where  the  parameters  are  as  defined  above. 

7.  S*,  the  source-group  Index,  Is  tested;  if  S*  Is  less 

•  n* 

than  (from  Library  3)»  (for  a  =  1  through  a^ax+l) 

Is  replaced  by  the  current  value  of  plus  and  h”* 

Is  replaced  by  the  current  value  of  h”*  plus  h“j  then  S* 

Is  Incremented  by  one,  and  the  routine  goes  to  Step  2. 

If  S*  =  3^ax«  current  value  of  h"*  Is  the  heat  genera¬ 
tion  rate  due  to  all  source  points,  the  current  values  of 
DjJ*  (a  r  1  through  a^ax^  spectral  points  of  the 

energy  spectrum  of  the  neutron  number  flux  due  to  all  the 

n« 

source  points,  and  Pamax4i  total  dose  rate. 

The  data  generated  by  this  routine  for  program  output  are: 

1.  Da  and  D^*  (for  a  =  1  through  a^ay )  -  the  spectra  of  the 
neutron  number  flux  summed  over  each  source  group  and 


over  all  source  groups,  respectively. 


2. 


^n*♦ 


"ilmax-M'  "Cix*!  - 


point,  per  source  group,  and  total,  respectively. 

n##  n  n* 

3>  H  ,  H  ,  H  -  the  heat  generation  rate  per  source  point, 
per  source  group,  and  total,  respectively. 
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Ill  UTILIZATION 


3.1  General 

This  section  lists  the  required  formats  for  library  and 
problem  Input  for  both  programs.  The  data  and  program-control 
parameters  required  for  utilization  of  these  programs  are  grouped 
Into  three  libraries  and  a  problem  deck.  In  addition*  there  are 
two  libraries  peculiar  to  the  simple-geometry  program  (C-17)«  one 
of  which  controls  readout  of  Information  from  the  library  tape* 
while  the  other  allows  deletion  of  decks  from  the  library  tape. 

The  Input  data  for  the  program  are  Input  In  blocks  called 
data  files.  These  data  files  contain  the  required  computational 
and/or  program-control  data  and  an  end-of-flle  symbol  which  consists 
of  an  asterisk  (*)  and  an  Identification  number.  The  first  piece 
of  data  In  a  file  must  start  on  a  new  card  and  any  alphabetlcs  In 
the  data  file  must  be  restricted  to  the  columns  shown  In  the  for¬ 
mats  (Figs.  3-1  through  3-8).  The  data  Is  restricted  to  Columns 
2  through  62  of  the  cards*  and  the  end-of-flle  symbol  should  begin 
In  Column  38.  (This  symbol  may  be  put  on  a  separate  card  If  de¬ 
sired.) 

Except  for  the  alphabetlcs*  the  conversion  digits  (the  6,  7 » 
or  8  which  must  be  put  In  Column  1)*  and  the  end-of-flle  symbol* 
the  data  are  not  restricted  to  any  particular  columns  within  the 
data  field;  rather*  only  the  order  of  the  Input  Is  Important*  and 
more  or  less  data  can  be  put  on  a  card*  If  desired*  than  Is  shown 
In  the  formats. 
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The  conversion  digit  7  Is  restricted  to  use  as  shown  In  the 
formats^  and  all  numbers  on  a  7  conversion  card  must  be  preceded 
by  a  sign  and  must  not  have  a  decimal  point;  this  conversion  Is 
used  for  control  parameters.  The  6  and  8  conversions  are  used  for 
data  Input  and  may  be  used  Interchangeably  from  card  to  card. 

For  6,  or  fixed-point  conversion,  numbers  are  Input  In 
normal  fashion  with  signs  and  decimal  points.  Signs  must  be  In¬ 
put  In  all  cases  and  occupy  a  column  In  front  of  their  respective 
numbers.  The  decimal  point  Is  assumed  to  be  to  the  right  of  the 
last  digit  If  It  Is  omitted.  A  fixed-point  number  may  consist  of 
up  to  11  digits,  with  never  more  than  10  digits  to  the  right  of 
the  decimal  point. 

For  8,  or  floating-point  conversion,  numbers  are  represented 
by  the  sign  of  the  number,  followed  by  from  one  to  eleven  digits 
representing  the  fraction,  followed  by  the  sign  of  the  exponent 
and  one  or  two  digits  representing  the  exponent.  Thus,  the  Input 
number  x  Is  defined  to  be 

X  s  F  •  10®  , 

where  O.lSF^l,  orF  =  0 

and  -  372  E  S38. 

The  above  restriction  on  E  Is  approximate,  since  this  restriction 
Is,  In  general,  dependent  on  the  number  of  digits  of  the  fraction 
F.  The  actual  restrictions  on  Input  numbers  az*et 
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1.  For  any  number  x,  |x  |Sl.70l4ll82l6  x 

2.  The  smallest  number  the  program  will  accept  Is 
shown  In  the  table  below  In  terms  of  the  number 
of  digits  In  the  fractional  portions 


Number  of 
Digits 

In  Fraction 

Minimum  Number 

Number  of 
Digits 

In  Fraction 

f— — — — 

Minimum 

Number 

Fraction 

Exponent 

Fraction 

Exponent 

1 

-37 

7 

t 1000000 

-31 

2 

-10 

-36 

8 

i 10000000 

-30 

3 

^100 

-35 

9 

-100000000 

-29 

4 

-1000 

-34 

10 

t 1000000000 

-28 

5 

-10000 

-33 

11 

^1000000000 

-27 

6 

liOOOOO 

_ _ 

Thus,  from  the  table  It  Is  seen  that  -^l  -  37  (  1.0  x  10”^®)  and 
-22  -  36  (-2.2x10"^^)  are  allowable  Input  numbers,  while  +-12  -  37 
(■♦•1.2x10"^®)  Is  not.  Note:  for  floating  point  conversion,  the 


decimal  point  Is  assumed  to  be  at  the  extreme  left  of  the  fraction, 
but  decimals  do  not  appear  In  this  representation. 

The  error  printout  from  the  programs  are  self  explanatory, 
but  the  errors  found  by  the  executive  program  (CF-109),  are  one 
of  the  following: 

Error  Number  1  -  The  program  has  found  an  undefined 

symbol.  Illegal  double  punchy  or 
Illegal  blank  column  In  the  data 
field  of  some  library  or  data  card. 

Error  Number  2  -  Exponent  trouble  In  floating  point. 

Error  Number  3  -  Undefined  control  punch  In  Column  1 

of  some  library  or  data  card. 


A  detailed  descrdptlon  of  the  library  and  data  formats  Is 
given  In  the  following  sections. 

Note:  The  Library  1  and  Library  2  fonnats,  and  hence  the 
data  for  these  libraries,  are  the  same  for  both 
programs . 

Note:  All  libraries  require  7  conversion  for  the  first 

data  file  in  the  library,  and  6  or  8  conversion 

for  all  other  data  files. 

3.2  Library  1:  Material  Data  (Fig.  3-1) 

This  library  is  composed  of  an  arbitrary  number  of  decks, 

each  deck  containing  all  of  the  data  pertaining  to  one  material. 

Two  types  of  decks  may  be  written  into  this  library.  The  first 

of  these,  which  Is  the  type  used  for  neutron  reference  materials, 

must  contain  a  *15  data  file  (neutron  differential  number- 

spectra  curve-fit  coefficients).  For  this  type  deck,  a  ^  0, 

max 

and  EG  9^0,  Each  of  the  other  data  files  of  this  library 

(Fig.  3-1)  may  or  may  not  be  listed.  The  second  type  of  deck 

is  that  used  for  non-reference  materials;  this  deck  must  contain 

a  *13  data  file  (neutron  heating  coefficients).  In  this  case, 

a^ax  =  0*  eg  s  0  and  no  *14  or  *15  data  files  may  be  put  in  the 

deck;  each  of  the  other  data  files  of  this  library  may  or  may 

not  be  listed,  A  non-reference  material  deck  may  be  made  Into 

a  reference  material  deck  by  an  appropriate  change  In  a  and 

nisix 

EG,  and  by  adding  a  data  file  *15.  Appendix  D  gives  typical 
Library  1  deck  lists  for  various  materials. 
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Columns 


NPC  I4.S24 
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'aMAX+1,1  ^AMAX+1,2  •  *  •  ‘^AMAX4-1,5 


Each  Library  1  deck  may  contain  the  following  data  files t 
A.  Control  Parameters  (*10).  This  data  file  consists  of  one 

card  (Plg. 3-1)  containing  the  following? 

1.  Library  type  Identification,  l.e,,  LIBRARY  DATA-#-l, 
(Columns  2  through  16) 

2.  M,  the  material  number  (a  positive  Integer).  This 
number  serves  as  a  deck  Identification  number,  and 
hence  the  decks  should  be  numbered  consecutively. 

This  number  Is  also  used  by  the  program  to  determine 
which  material  deck  Is  to  be  Identified  with  a  given 
region  for  a  particular  problem  (Pigs.  3-7  and  3“8)* 

3.  EO,  neutron-energy -mode  Identifier.  This  number  Is 
Input  as  an  Integer,  0£E0£3>  Each  non-zero  value  of 
EQ  Identifies  a  particular  set  of  neutron  energies, 
the  first  four  of  which  are  listed  In  Table  3-1 •  The 
number  EO  Is  non-zero  only  If  the  material  deck  contains 
a  *15  data  file  (neutron  differential  number-spectra 
curve-fit  coefficients)  and.  In  this  case,  EQ  refers 

to  the  set  of  neutron  energies  (from  Table  3-1)  to  which 
the  neutron  curve  fits  correspond. 

4.  %ax«  number  of  energies  for  which  the  neutron 

differential  number  spectrum  is  defined  In  the 
material  deck,  non-zero  for  reference  materials 

only  and.  In  general,  0  f  -10. 
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Neutron  Energy  Intervals  and  Histogram  Factors 
for  the  First  Four  Neutron-Energy  Modes 
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B.  Removal  Croas  Sections  (^11) 

Past  neutron  removal  cross  sections  are  In  units  of  cm^/gm. 
(The  cross  sections  for  the  various  material  decks  are  from 
Reference  8.) 

C .  Qamma  Data  (*12) 

If  a  Library  1  deck  contains  a  *12  data  file,  then  each  of 
the  elements  listed  below  must  be  put  In,  even  though  these 
data  are  Input  as  zeros.  These  data  are  Input  In  order  of 
descending  energy  and  the  energy  sets  used  are  listed  In 
Table  3*2.  The  coefficients  In  this  data  file  are* 

1.  through gamma-ray  attenuation  coefficients 
or  mass  absorption  coefficients  In  units  of  cm^/gm. 

These  coefficients  have  been  Interpolated  trcm 
Reference  6. 

2.  Zeff^  through  Zeff^^ij,  effective  atomic  number.  For 
elemental  materials,  these  are  the  atomic  number  of  the 
material,  hence,  are  Independent  of  energy;  but  for 
mixtures  and  compounds,  Zeff  Is  found  by  first  computing 
the  absorption  coefficient  per  electron  (using  the 
foznnulae  below),  then  from  the  computedyA^  an  Interpola¬ 
tion  of  Z  Is  made  as  a  function  of  the  absorption 
coefficient  per  electron  for  the  elements.  This  Zeff 
Is,  In  general,  a  function  of  photon  energy. 

The  absorption  coefficient  per  electron  for  mixtures 
and  compounds  Is  given  by 
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N 


Is  the  number  of  elements  in  the  material. 


1 

—W -  ^ 

2: 

1=1 

the  weight  fraction  of  the  1^  element  for  mixtures 

niA^ 

—  where  n^^  equals  the  number  of  atoms  of  the 

^  ^  1^  kind  (which  has  atomic  weight  A^) 

in  a  molecule  for  compounds , 

2 

is  the  absorption  coefficient  (cm  /gm) ,  and 


-Z 
—  i 


is  the  number  of  electrons  per  gram  of  the  1^  material. 


3. 


^  through  H,-,  ,  gamma-ray  flux-to-heat  conversion 


X 


factors . 
ZH 

CV  1  1 


These  factors  are  given  by 
e  si^  ^  (1.602  X  10"^3  watts/Mev), 


where 


^1 


the  photon  energy  for  which  the  coefficients 
are  being  computed  (in  Mev), 

is  the  material  electron  density  (in  electrons/gm), 
is  the  Kleln-Nlshlna  scattering  cross  section 
as  stated  on  page  25  of  Reference  9  (in  cm^/ 
electron)  in  hydrogen  for  the  energy  E^,  and 
is  the  gamma-ray  attenuation  coefficient  for 
the  1^  energy,  as  defined  above. 


D,  Neutron  Plux-to-Heat  Conversion  Factors  (*13) 

As  stated  above,  If  the  Library  1  deck  pertains  to  a  non¬ 
reference  material,  the  deck  must  contain  a  *13  data  file 
(this  may  be  simulated  by  a  card  containing  one  or  two  zeros 
and  the  end-of-flle  symbol,  *13).  This  data  file  consists 
of  from  one  to  five  sets  of  fast-neutron  flux-to-heat  con¬ 
version  factors  and  an  Identifying  symbol  as  defined 

below. 

1.  ®hnaxEOi  number  of  energies  for  which  the  neutron 

differential  number  spectrum  Is  defined  for 
neutron  energy  mode  1. 


®maxEQ2  ®re  similarly  defined  for  neutron 

energy  modes  2  through  5* 

r  -  M  « 

2.  through Ort  fast-neutron,  flux-to-heat  con- 


CL 


'■mmx 


version  factors.  These  factors  are  proportional 
to  group-averaged  neutron  elastic-scattering 


cross  sections  (or  other  applicable  fast-neutron 
cross  sections),  and  are  given  by 


where 


^%*1  neutron  scattering  cross  section  per  unit 

density  averaged  over  energy  group  of  the 

neutron  energy  mode  1,  and 


^ is  the  fraction  of  the  Initial  neutron's  energy 
dissipated  as  heat. 
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A  more  thorough  treatment  of  these  coefficients,  and  the 
coefficients  for  other  neutron  reactions  Is  given  In  Appendix  B 


3  •  J  ^ -  *  • '  hrough  6 

*'*  . 
past-neutron  flux-to-heat  conversion  factors  for  neutron 

energy  modes  2  through  3$  respectively.  The  units  of 


will  be  wat<t  cm^/neutron-gm,  l.e.,  the  heat  In  watts^ 
generated  by  one  neutron  in  traveling  one  centimeter  in  a 
material  of  unit  density. 

E.  Neutron  Reference-Material  Data  (*14  and  *15) 

The  last  two  data  files  (*14  and  *15)  shown  in  Figure  3-1  must 
be  listed  for  all  reference  materials. 

The  first  of  these  data  files  (*14)  consists  of  the  neutron 
histogram  factors;  these  are  the  numerical  integration  factors 
used  In  the  Integration  over  neutron  energy  required  for  the 
neutron  heating  calculation.  This  data  file  consists  ofi 


hf  ulth rough  .hf  ,  neutron  hlstograni  factors.  These 

^max 

factors  are  listed  only  for  reference  materials,  and  the 
energy  mode  used  to  derive  these  factors  Is  that  mode  used  to 
define  the  neutron  differential  number  spectra  for  the  material 
for  which  the  deck  Is  defined. 

The  second  data  file  (*15)  contains  the  curve-fit  coefficients 
for  neutron  differential  number  spectra.  In  general.  It  was 
not  possible  to  fit  these  spectra  to  one  set  of  coefficients 
and,  hence,  two  sets  are  required  along  with  a  break  point 
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BP  (If  the  equivalent  thickness  t._5BP«  the  first  set  Is  used, 

eq 

and  If  t^q^BPf  the  second  set  Is  used).  This  data  file  con¬ 
sists  oft 

1.  C,  ,  through  C,  _  the  set  of  curve-fit  coefficients  for 

the  neutron  differential  number-spectra  point  for  the 

highest  energy.  Two  sets  of  these  coefficients  are  needed 

for  each  energy;  one  set  for  0  -  t  -  BP«  the  second  for 

eq 

BP  —  t_„  ^  250  gpn/cm^.  The  coefficients  are  for  the 
eq 

function 

P(t)  ;  exp(C]^t^-f  Cgt^  4  ‘h  C4t  +  C^) 

where  t  Is  to  be  In  gm/cm^.  ^2,1  ^2,3,  ®3,1 

through  Co  c  ...«  C.  through  C.  are  the  sets  of 

nnax,l  “max,  5 

curve-fit  coefficients  for  the  remainder  of  the  neutron 
differential  number  spectra  and  are  In  order  of  decreas¬ 
ing  energy.  The  set  of  coefficients  C.  through 

“maxfl, 1 

C  fit  the  fast-neutron  dose  rate  In  the  same 

®max'fl,5 

manner  as  for  the  neutron  spectra. 


Note:  All  of  the  neutron  dose-rate  curves  were  computed 
using  the  flux-to-dose  conversion  factors  set  forth  In 
NYO-6269  (Ref.  10).  The  curve-fit  coefficients  must  be 
Input  such  that  che  Independent  variable  t  will  have 
units  of  gm/cm|  and  the  function  F(t)  will  have  units 
of  neutrons/cm“-sec-Mev  per  Incident  neutron  for  the 
spectral  curves  and  mllllrem/hr  per  Incident  neutron 
for  the  dose-rate  curves.  The  data  from  which  these 
coefficients  were  generated  are  given  In  References  11,  12, 
and  13. 
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2.  BP,  curve-fit  break  point.  This  number  is  used  to  indi¬ 
cate  which  one  of  the  two  sets  of  curve-fit  coefficients 


for  the  neutron  spectra  are  to  be  used  to  generate  the 
spectra  for  a  particular  value  of  t,  since  one  set  of 
coefficients  is  defined  for  t^BP,.  and  the  second  for 
tiBP.  Note  that  BP  has  the  same  units  as  t;  namely, 
gm/cm^. 

3.3  Library  2:  Gamma  Data  (Fig.  3-2) 

This  library  is  composed  of  one  deck  which  contains  gamma- 
energy  values,  flux-to-dose  conversion  factors,  coefficients  for 
the  curve  fits  to  the  differential  energy  spectra,  histogram  factors 
required  for  both  the  gamma  spectral  Integration  and  the  integra¬ 
tion  over  energy  required  in  the  gamma  heating  and  dose-rate 
calculations,  and  the  coefficients  for  the  curve  fits  to  the  edge 
corrections  of  Berger  and  Doggett  (Ref.  7). 

Appendix  D  gives  the  Library  list  presently  in  use  at  GD/PW. 

Library  2  is  required  to  have  a  data  file  *20  (edge  correc¬ 
tions);  and  may  contain  the  following  data  files: 

A.  Control  parameters  (*16) 

Library-type  Identification;  i.e.,  LIBRARY  DATA  h2 
(Columns  2  thru  I6) 

B.  Gamma-ray  initial/final  energies,  and  flux-to-dose 
conversion  factors  (*17) 

through  £^14*  degraded  gamma  energies.  The  four¬ 
teen  values  are  listed  in  Table  3*2. 

2.  through  £514,  initial  gamma  energies.  These  are 

the  same  as  the  degraded  energies  with  the  exception 
of  the  lowest  energy  (0.25  Mev),  which  is  deleted. 
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_2(D  ^YD 

3.  6.  through  <5  ,  flux-to-dose  conversion  factors, 

1  ft,  hiO'X 

These  are  listed  in  Appendix  D. 

Qamma-ray  differential  energy  spectra  coefficients  (*l8) 

1.  Zeff(maxi),  break  point  between  set 

2.  through  ''104,6' 

coefficients  for  the  curve  fits  to  the  gamma-ray 
differential  energy  spectra  for  the  case  0  SZ—  26, 
These  coefficients  are  for  the  function 


§ 


where 

is  the  gamma  absorption  coefficient 
at  energy  and 

f(T)  is  the  gamma-ray  differential  energy 
spectra  as  defined  in  Reference  2, 

The  dependence  of  the  A's  (as  shown  by  the  second 

subscript)  on  the  energy  Indices  a  and  b  (or 

equivalently  on  degraded  and  initial  energy)  is  as 


shown  in  Table  3-2, 

3.  Zeff(max2),  break  point  between  set  Ag  and  set 

A-  „„  of  the  A's. 
xx 

4.  The  sets  of  coefficients  A?  .  through  A?  , 

1|1  1^0 

*104,1  *104,6  *1,1  *^04,1 

3 

^1Q4,6  similarly  defined  for  the  ranges 
26  —  Z  ~74  and  7^—  92,  respectively. 


5.  Zeff(max3)j  a  value  as  large  or  larger  than  any  Zeff 
to  be  used. 

Notet  There  are  three  *18  files,  one  for  each  set 
Zeff (max)  and  A's. 

D.  Histogram  Factors  (*19) 

hfy^  through  histogram  factors  for  the  spectral 

Integrations  and  for  the  Integration  over  degraded  photon 
energy  required  for  the  gamma  heat  generation  rate  and 
dose-rate  computations.  These  factors  are  Input  In  the 
same  order  (as  regards  Initial  and  final  energy  Indices 
data)  as  the  differential  energy  spectra  (*18)  file. 

E.  Edge  Corrections  (*20) 

1.  ^r(maxl),  breakpoint  between  set  B,  _  and  _ 
of  the  B's. 

2.  through  through  tho 

coefficients  for  the  curve  fits  to  the  edge 
corrections  to  the  gamma-ray  differential  energy 
spectra  for  the  case  0  5>trl.  4.  These  coefficients 
are  for  the  function 

where 

b  is  the  gamma-ray  absorption  coefficient 
at  energy  E|^  , 
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g(T)  Is  edge  correction  of  Reference  7>  l.e., 
g(T)  - 

with  the  buildup  factors  in  an  infinite 
medium,  and  for  a  slab,  B(t,t), 

as  defined  in  Reference  4. 

3.  through 

similarly  defined  for  4'^^r  <20.  The  first 
subscript  on  the  B's  refers  to  initial  photon 
energy  and  is  the  energy  counter  b,  as  used  in 
(C)  above. 

4.  yVr(max2),  a  value  larger  than  any  to  be  used. 
3.4  Library  3t  Source  Data 

This  library  is  composed  of  an  arbitrary  number  of  decks, 
each  deck  containing  the  source  data  for  a  specific  problem. 

Each  deck  contains  from  one  to  six  source  groups  and  the  required 
deck-control  data.  Each  source  group  contains  the  source-group 
gamma  spectrum  (13  quantities  which  define  the  gamma  spectrum  at 
each  source  point  in  the  group)  and  four  quantities  for  each 
source  point  in  the  group,  namely,  3  source-point  coordinates 
and  the  source  Intensity.  Each  deck  is  limited  to,  at  most,  6 
source  groups  and  to  a  total  of  500  source  points  which  may  be 
divided  among  the  source  groups  in  any  manner.  This  library  is 
required  to  have  a  data  file  *22. 
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3.4.1  Library  3  For  Simple- Qeome try  Program  (Fig.  3-3). 

Library  3  for  the  simple -geometry  program  (C-17)  con¬ 
tains  the  following  data  files: 

A.  Control  parameters  (*21) 

1.  Library-type  Identification,  l.e.,  LIBRARY  DATA-#-  3 
(columns  2  through  l6) 

2.  N  ,  the  number  of  source  groups,  li  N.  —  6. 

c  c 

3.  ^s2*  ***'  ^sNc'  oumber  of  source  points  In  the 

1st,  2nd,  source  group,  respectively. 

B.  Source  data  (*22).  There  must  be  N  of  these  fields, 

c 

one  for  each  source  group. 

1.  Sj,  Sg#  ...»  ^13'  gamma-ray  spectrum,  tabulated 
In  order  of  decreasing  Initial  photon  energy  (from 
10  Mev  to  0.5  Mev). 

2.  X,  y,  z,  I,  the  x,  y,  and  z  coordinates  of  the 
source  point  and  the  source  Intensity;  one  set 

of  these  four  numbers  must  be  Input  for  each  source 
point  In  the  source  group. 

Note:  For  reactor  radiation  problems,  or  any  other 
problem  where  the  neutron  portion  of  the  pro¬ 
gram  Is  to  be  used,  the  gamma  spectra  should 
be  normalized  to  the  number  of  photons  per 
fission -Mev  so  that  the  source  Intensity 
may  have  units  of  fissions  per  watt  of  reactor 
power.  If  the  program  Is  to  be  used  for 
secondary -gamma  calculations,  or  calculations 
of  a  similar  nature  where  the  neutron  portion 
of  the  program  Is  not  used,  the  normalization 
of  the  spectra  Is  arbitrary,  since  the  gamma 
spectra  multiplied  by  the  source  Intensity 
must  result  In  units  of  photons/sec -Mev  for  a 
given  power  level. 
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NPC  I4.S2« 
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3.4.2  Library  3  for  Comp lex- Geometry  Program  (Fig  3*4) 

A  library  3  for  the  complex- geometry  program  (L-63)  con¬ 
tains  the  following  data  files  s 

A.  Control  parameters  (*21) 

1.  Library-type  Identification,  l.e.,  LIBRARY  DATA  -4-^  3 
(columns  2  through  l6). 

2.  N  ,  (defined  In  Section  2.4.1,  A) 

c 

3.  Nai,  Ng2*  ...>  Ngjjjj,  (defined  In  2.4.1,  A) 

4.  CT2#  ...»  CTjj,  coordinate  type 

If  CT  ;  1,  the  source  points  are  to  be  Input  In 
cartesian  coordinates. 

If  CT  ;  2,  the  source  points  are  to  be  Input  In 
cyllndrlc  coordinates. 

5.  CS, ,  CS^,  ...,  CS„  ,  coordinate  system  Identification 

12  Nc 

number.  Each  of  these  numbers  must  be  zero  (base 
coordinate  system)  or  must  correspond  to  the  COORN 
In  the  problem  deck  which  defines  the  coordinate 
system  In  which  the  source  points  for  the  particular 
group  was  defined. 

B.  Source  Data  (*22).  There  must  be  N  of  these  fields,  one 
for  each  source  group. 

1 .  ^  ^2^  ^3^  •••»  ^13^  (defined  In  3.4.1,  B) 

2.  X,  Y  or  0,  Z  or  R,  I,  the  x  coordinate,  and  either 
the  Y  and  Z  coordinates  (if  CT  =  1  for  this  source 
group),  or  the  6  and  R  coordinates  (If  CT  r  2  for 
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LIBHARY  DATA  (LTYPE  =  3) 
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this  source  group),  and  the  source  Intensity.  One 
set  of  four  numbers  must  be  Input  for  each  source 
point  in  the  source  group.  All  source  points  in  the 
source  group  must  be  in  either  Cartesian  or  cylindric 
coordinates,  but  not  both. 

The  note  above  for  the  simple  geometry  program  also 
applies  here,  giving  the  relationship  between  the 
^s  and  I *8. 

3.5  Deletion  and  Llstout  Libraries  for  the  Simple -Geometry  Program 

The  two  libraries  described  below  are  used  with  C-17  only. 

3.5.1  Library  4|  Library  Deletion  List  (Pig.  3-5) 

The  use  of  this  library  is  optional,  and  it  is  used 
only  to  list  those  decks  which  are  to  be  deleted  from  the  libraries. 
The  numbers  listed  in  this  library  are  those  which  appear  In 
Columns  63  through  69  of  the  decks  to  be  deleted.  (The  numbers 
must  have  plus  signs  associated  with  them).  The  alphabetic 
DELETE  must  appear  in  Columns  2  through  7  of  the  first  cai?d  of 
this  library.  Fifty  or  less  decks  may  be  deleted  by  using  one 
Library  4. 

3.5.2  Library  5:  Deck  Listout  (Fig,  3-6) 

The  use  of  this  library  is  optional,  and  it  is  used 
only  to  list  those  decks  which  are  to  be  printed  out  from  the 
library  tapes.  The  numbers  listed  in  this  library  are  those 
which  appear  in  Columns  63  through  69  of  the  decks  to  be  printed 
out  (the  numbers  must  have  plus  signs  associated  with  them). 
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The  alphabetic  LPRINT  must  appear  In  Columns  2-7  of  the  first  card 
of  this  library.  Fifty  or  less  decks  may  be  listed  using  one 
Library  5* 

3.6  Card  Identification  Field  for  Library  Decks 

1.  Columns  63  through  66  are  left  blank.  (These  are  to 
contain  the  Job  number  specified  by  the  Computing 
Laboratory) . 

2.  Column  67  contains  the  last  digit  in  the  year,  l.e,, 
0  for  i960,  1  for  1961,  etc. 

3.  Columns  68  and  69  contain  the  deck  number;  starts 
with  01  for  any  given  set  of  library  Input. 

4.  Columns  70  through  72  contain  card  numbers;  starts 
with  001  In  any  given  deck. 

5.  Column  73  contains  an  L. 

6.  Columns  78  through  80  contain  C17  or  li63> 

3.7  Problem  Deck 

This  deck  contains  the  program  options  and  that  data  peculiar 
to  the  problem  at  hand,  namely,  the  parameters  required  to  define 
the  problem  geometry,  the  detector  coordinates,  the  material 
numbers  of  the  neutron  reference  materials  to  be  used,  and,  for 
L-63,  the  coordinate  systems  used  to  define  the  geometry,  source 
and  detector  points. 

The  geometry  for  these  programs  Is  defined  by  using  a  con¬ 
sistent  set  of  x-planes  to  represent  each  region  In  the  problem 
after  the  method  set  forth  In  Section  2.1.  In  this  context,  for 
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the  simple  geometry  program,  an  x-plane  is  taken  to  be  an  x-value 
and,  for  the  Cartesian-geometry  option,  two  Y  and  two  Z  values 
(these  five  numbers  are  sufficient  to  determine  a  particular 
rectangle  with  respect  to  the  assumed  coordinate  system),  or,  for 
the  cylindrlc  geometry  option,  two  R  values  (these  3  numbers  are 
sufficient  to  determine  a  particular  circular  annulus  with  respect 
to  the  assumed  coordinate  system).  Thus,  two  adjacent  x-planes, 
assuming  linear  interpolation-  or  circular  interpolation  for  the 
spherical  option  in  L-63,  between  analogous  Y  and  Z  values  or 
R  values...  form  a  volume  element,  and  each  region  is  built  up  of 
these  volume  elements.  Similar  definitions  hold  for  the  complex- 
geometry  program. 

3.7.1  Problem  Deck  for  the  Simple -Geometry  Program  (Fig.  3-7) 
The  problem  deck  for  the  simple-geometry  program  (C-17) 
contains  the  following  data  fields  1 

Note:  In  this  deck  the  *6  and  *8  data  files  require 

7  conversion;  all  other  files  use  either  6  or 

8  conversion. 

A.  Control  parameters  (*6) 

This  file  contains  the  following: 

1.  Deck  identification,  l.e.,  the  alphabetlcs  PROBLEM 
"  DATA  must  appear  in  Columns  2  through  13. 

2.  IDg  and  ID^,  the  identification  of  the  decks  of 
Libraries  2  and  3  to  be  used  in  this  problem.  These 
are  the  numbers  contained  in  Columns  63  through  69  of 
these  libraries . 
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Card  Columns 


3.  the  number  of  geometric  regions  defined  for  the 
problem.  This  Includes  the  dummy  region  which  must 
always  be  defined  and  must  be  the  last  region. 

4.  the  number  of  neutron  calculations  to  be  made 

or,  equivalently,  the  number  of  neutron  reference 
materials  listed  (by  material  number  M)  In  the 
problem  deck.  r  0  If  no  neutron  calculation  Is 

to  be  made,  and.  In  general,  0  i  i  10. 

3.  U,  gamma  calculation  option. 

U  ;  1  -  differential  energy  spectra  for  an 
Infinite  media  used 

U  r  2  -  edge -corrected  differential  energy  spectra 
used 

U  2  3  "  both  calculations  (1  and  2)  made. 

Volume  Parameters  (*7) 

This  field  contains  all  of  the  parameters  required  to 
define  one  volume  of  the  geometry.  Each  of  these  fields, 
excepting  the  last  one  (which  defines  the  dummy  region), 
must  contain  the  following: 

1.  Region  Identification.  The  alphabetlcs  REON  must 
appear  In  Columns  2  through  5  of  the  first  card  In 
a  *7  field,  and  Columns  6  and  7  of  this  card  must 
contain  a  two-digit  region  identification  number 
(the  Nj^  regions  must  be  numbered  In  sequential 


order  as  they  appear  In  the  deck,  l.e.,  01,  02, 


2.  M«  material  number.  The  number  of  the  Library  1 
deck  which  corresponds  to  the  material  for  which  the 
region  Is  defined. 

3.  ^  $  the  density  of  the  material  M  In  the  region 

(gm/cm^). 

4.  the  number  of  x-planes  used  to  define  the  region. 

5.  K,  geometry-type  parameter. 

K  =  1  -  region  Is  cyllndrlc 
K  2  2  -  region  Is  Cartesian 

6.  uncertainty  parameter  (cm).  The 

stepplng-polnt  method  precludes  an  exact  determination 
of  the  points  at  which  the  "source-detector  line" 
Intersects  the  region  bounds.  Hence«  It  Is  necessary 
to  Input  a  number  to  define  the  greatest  allowable 
error  In  the  boundary  determination.  Since  the 
stepplng-polnt  method  Is  an  Iterative  technique^  an 
excessively  small  value  of  will  require  a  large 

number  of  Iterations  and,  hence,  will  Increase  the 
machine  time  required  to  run  the  problem.  On  the 
other  hand,  too  large  a  value  for  will  result 

In  excessive  error  In  the  boundary  determination 
and,  hence,  the  final  results  will  be  In  error  due 
to  these  uncertainties  In  the  determination  of  the 
boundaries . 
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One  method  of  determining  K  ^  is  to  require  that  it 

min 

be  smaller  than  some  fraction,  say  0,01,  of  the 
smallest  relaxation  length  (either  neutron  or  photon) 
of  the  particles  in  the  materials  of  the  system  being 
considered.  For  example,  consider  a  system  composed 
of  the  GTP  and  a  shield  of  water,  aluminum.  Iron  and 
Portland  concrete.  Prom  Reference  8,  the  removal 
cross  sections  for  the  above  are:  0.09^9  cm“^, 

0.101  cm^^,  0.0788  cm“^,  O.I688  cm"^.  and  0.080I  cm“^, 
respectively.  Prom  Reference  14,  the  largest  values 
of  the  photon  attenuation  coefficients  (In  each  case 
for  a  photon  energy  of  O.25  Mev)  are:  0.225  cm"^, 
0.126  cm“^,  0.302  cm“^,  0.918  cm"^,  and  0,270  cm~^, 

respectively.  The  largest  of  these  numbers  (the 

attenuation  coefficient  for  E  ;  O.25  Mev  In  Iron) 

Is  0.919  cm~^  which  corresponds  to  a  relaxation 

length  of  1.089  cm.  Hence,  If  K  .  -  .01  x  I.O89  - 

min  •  “* 

0.01089  cm,  then  the  boundaries  will  be  determined 
so  that  the  error  Involved  In  determining  the  source- 
detector  path  length  In  any  region  will  be  less  than 
0.01  of  a  relaxation  length  for  either  photons  or 
neutrons  In  this  region. 

7.  The  following  data  (in  cm)  must  be  entered  for  each 
of  the  x-planes  used  to  define  the  region: 
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a.  X  $  the  x-plane  x-ooordlnate 

p 

Note:  the  x-planee  for  a  given  region  must  be 
listed  In  order  of  Increasing  Xp. 

b.  Either  and  Z^.  the  i-plane 

bounds  for  a  cartesian  region  (If  K  =  2),  or 
and  R^,  the  x-plane  bounds  for  a  cyllndrlc 
region  (If  K  =  1). 

The  last  region  (or  dummy  region)  contains  only  Items  1«  2, 
and  3  above. 

C.  Neutron  reference  material  numbers  (*8) 

This  field  contains  material  numbers,  each  of  which 
corresponds  to  a  material  deck  containing  a  *13  data 
field.  (This  field  must  contain  the  alphabetic  REPMAT 
In  Columns  2  through  7  of  the  first  card  of  the  field). 

D.  Detector  and  program  option  parameters  (*9) 

This  data  file  contains  the  following} 

1.  Field  Identification.  The  alphabetic  DETECT  must 
appear  In  Columns  2  through  7  of  the  first  card  of 
this  field. 

2.  Np,  the  number  of  detector  points  for  the  problem. 

3.  The  following  data  are  Input  for  each  of  the 
detector  points: 

a.  x^,  y^,  and  detector  point  coordinates 

b.  X,  the  detector  option 

X  r  1  -  Spectrum  and  dose  will  be  calculated. 

X  -  2  -  Spectrum,  heat,  and  dose  will  be 
calculated. 

X  2  3  Spectrum  and  heat  will  be  calculated. 
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c.  Y,  the  radlatlon-type  option 

Y  :  1  -  Geometry  print.  The  source-detector 

distance  and  penetration  distance 
through  all  non-void  materials  will 
be  printed  for  each  source-detector 
pair. 

Y  -  2  -  Gamma-ray  data  only  will  be  computed. 

Y  =  3  '  Gamma-ray  and  neutron  data  will  be 

computed. 

Y  3  4  -  Neutron  data  only  will  be  computed. 

d.  Q,  the  source  point  option 

Q  r  1  -  Data  for  each  source  point  as  well 

as  data  summed  over  each  source  group 
and  over  all  source  groups  will  be 
printed. 

Q  =  2  -  Data  summed  over  each  source  group 
and  over  all  source  groups  will  be 
printed . 

3.7.2  Problem  Deck  for  the  Complex- Geometry  Program  (Pig.  3-8) 
The  problem  deck  for  the  complex  geometry  program  (L-63) 
contains  the  following  data  file: 

A.  Control  parameters  (*6) 

This  field  contains  the  following: 

1.  Deck  identification  (as  defined  In  3»7.1-A) 

2.  IDg  and  ID^  (as  defined  In  3.7.1-A) 
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MASTKR  ^MASTN-NMR)  (Mdummy )  {  ^  dummy ) 


FIGURE  3-8.  PROBLEM  DECK  FOR  THE  COMPLEX -GEOMnRY  PROGRAM 


3.  the  number  of  master  regions  defined  for  the 
problem.  This  Includes  the  "dummy  volume"  which  must 
be  listed  as  a  master  region  and  also  must  be  the  last 
master  region. 

4.  Npj,  (as  defined  In  3.7.1-A) 

5.  U  (as  defined  In  3«7.1-A) 

6.  the  number  of  coordinate  systems  used  to 
define  the  problem  geometry  (the  base  system  Is  not 
counted) . 

Note*  This  data  file  requires  7  conversion. 

B.  Volume  parameters  (*7) 

This  field  contains  all  of  the  parameters  required  to  de¬ 
fine  one  volume  of  the  geometry.  Each  of  these  fields 
excepting  the  last  one  (the  dummy  volume  which  must  be 
listed  as  a  master  region)  must  contain  the  following: 

1.  The  alphabetic  Identification 

MASTER  for  master  regions « 

REGION  for  regions,  or 
SUBREG  for  subregions 

In  Columns  2  through  7  of  the  first  card  of  the  file. 

2.  N  ,  the  volume  number,  which  must  be  preceded  by 

R 

a  plus  sign.  Volumes  must  be  numbered  consecutively 
by  type;  that  Is,  the  first  master  region  Is 
numbered  1,  the  second  +2,  etc.  Further,  If  some 
master  region  contains  regions,  the  data  files  for 
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these  regions  must  follow  the  data  file  of  the 
master  region  In  which  they  are  contained,  and 
these  regions  are  numbered  starting  with  +1. 

If  some  region  contains  subregions,  the  same  pro¬ 
cedure  Is  followed,  that  is,  the  data  files  for 
the  subregions  follow  the  data  file  for  the  region 
and  they  are  numbered  consecutively,  starting  with 
+1. 

3.  M  (as  defined  In  3.7.1-B) 

4.  Njj  (as  defined  In  3.7.1-B) 

5.  CT,  the  coordinate  type  used  to  define  the  volume. 

This  parameter  Is  as  defined  in  Section  2.3.2,  or 

10  b  below.  A  recapitulation  of  the  values  for  CT  is: 
CT  =  0  Simple  Cartesian  cross  section, 

CT  =  1  Complex  Cartesian  cross  section, 

CT  =  2  Simple  cyllndrlc  cross  section, 

CT  z  3  Complex  cyllndrlc  cross  section, 

CT  r  4  Rectilinear  cyllndrlc  cross  section, 

CT  =  5  Combination  I  cyllndrlc  cross  section, 

CT  =  6  Combination  II  cyllndrlc  cross  section. 

6.  CS,  the  number  of  the  coordinate  system  used  to  define 
the  volume.  CS  must  correspond  to  one  of  the  COORN 

of  problem  data  file  *10. 

7.  NREQ  or  NSUB,  the  number  of  volumes  contained  In  the 
volume  being  considered.  If  the  volume  Is  a  master 

ir- 


region,  then  NREQ  is  the  number  of  regions  contained 
In  the  master  region;  the  subregions.  If  any,  con¬ 
tained  In  one  of  these  regions  are  not  counted.  If 
the  volume  Is  a  region,  then  NSDB  Is  the  number  of 
subregions  contained  In  the  volume.  In  any  case, 
O^iNREQ  —10,  and  0  —  NSUB  ^10.  If  the  volume  Is 
a  subregion,  this  Item  Is  omitted. 

Note:  The  above  Items  require  7  conversion  and, 
hence,  must  not  appear  on  the  same  card  as 
any  of  the  following  Items  (which  may  have 
either  6  or  8  conversion). 

8.  /^,  the  density  of  the  material  M  for  the  volume 
(gm/cm^) . 

9.  (as  defined  In  3.7.I-B) 

10.  The  following  data  must  be  entered  (in  cm)  for  each 
of  the  x-planes  used  to  define  the  volumes: 

a.  Xp,  the  x-plane  x-coordlnate 

Note:  the  x-planes  for  a  given  volume  must  be  listed 
In  order  of  Increasing  Xp. 

b.  Set  of  numbers  required  to  define  an  x-plane 
cross  section.  These  numbers,  for  each  case, 
are  listed  below  In  the  order  in  which  they  are 
to  appear  In  the  program.  Along  with  each  set 
of  numbers  there  Is  also  shown  a  sketch  of  a 
typical  cross  section  which  they  are  to  repre¬ 
sent,  the  required  value  of  CT,  and  the 
restrictions  on  these  numbers. 


113 


Note:  For  all  of  the  figures  below,  the  positive 

x.-dlrectlon  is  assumed  to  be  out  of  the  paper. 


Simple  cartesian  (CT  =  0) 


Input  values  for  each  x-plane: 

^mln*  ^max'  ^mJn*  ^max* 
Restrictions : 

The  defined  area  is  a  rectangle 
whose  sides  are  parallel  to 
either  the  y  or  z  axis; 

^min  ^  ^max  *  ^mln  ^  ^max 


Y  Y 

min  max 

> _ I 

^max 

^mln 

- Y 


Complex  Cartesian  (CT  =  1) 


Input  values  for  each  x-plane: 


Yl^Zi  ,  Y2,Z2  ,  YyZ^  , 


Y4,  Z4  . 

Restrictions: 

(On  the  points  =  (Xp,  Yj^,  Z^) 
for  1  r  1#  2,  3#  and  the  sides 
PlPj  which  is  the  side  with  end 
points  Pj^  and  Pj) 

Orientation:  Neither 


Z 


P1P2  nor  P2J’4  is  parallel  to 


the  z  axis 

Neither  P^^P^  nor  ^2^*3  parallel  to 
the  y  axis 


Connectivity:  P^  is  adjacent  to  P2  and  P^  in  the 
sense  shown  in  the  figure  above. 
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Convexity: 


All  Internal  angles  must  be  less 
than  180°. 


Note:  For  all  of  the  cyllndrlc  regions  defined  below, 
the  following  restrictions  apply: 

1.  The  cylinder  must  be  coaxial  with  the 

X  axis  of  the  coordinate  system  In  which 
it  Is  defined. 

2.  All  angles  are  to  be  measured  In  degrees 
and 

360° 


Simple  cylinder  (CT  z  2) 


Complex  cylinder  (CT  =  3) 
Input  values  for  each  x-plane: 

%n*  **o»  ®1'  ®2* 

Restriction: 

%n  ^  ^o* 
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If  CT  is  preceded  by  a  minus  sign,  the  program  will  generate 
portions  of  spheres  with  the  cross  sections  above. 

B'.  Volume  parameters  for  the  dummy  volume  (*7) 

The  last  volume  defined  for  a  problem  must  be  the 
dummy  volume,  the  data  field  for  this  volume  must 
contain 

1.  Alphabetic  Identification,  MASTER,  in  Columns  2 
through  7 

2.  The  volume  number  which  must  be  (see  A, 3  above) 

3.  M,  (defined  in  3.7.1-B) 

4.  (defined  in  B^8  above) 

The  parameters  above  are  all  that  is  required  to  define  the 
dummy  volume,  but  it  should  be  noted  that,  while  7  conversion  is 
required  for  the  first  card  of  other  volume  data  file,  6  or  8  con¬ 
version  is  required  for  the  dummy  volume. 

C.  Neutron  reference  material  numbers  (*8) 

This  data  file  is  the  same  as  the  one  for  the  simple- 
geometry  program  (3.7.1-C). 

D.  Detector  and  program  option  parameters  (*9) 

This  data  file  contains  the  following: 

1.  Alphabetic  identification,  DETECT,  in  Columns  2 
through  7  of  the  first  card. 

2.  d*jjj^the  number  of  detector  groups 

3.  Dj,  Dg,  ...,  the  number  of  detector  points  in 

each  detector  group. 
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4,  CT^,  CTp,  CT^  ,  the  type  of  coordinate  used  to 

1  C. 

define  the  detectors  of  the  various  groups,  l.e., 

CT  =  1,  Cartesian  coordinates 
CT  =  2,  cyllndrlc  coordinates 

5.  CSj,,  CSjj  ,  the  coordinate  system  used  to 

tn 

define  the  detectors  of  the  various  groups;  each 
CS  ^0  must  correspond  to  some  COORN  of  data  file 
*10  below. 

Note:  The  above  listed  parameters  compose  the  first  *9 
data  file  in  a  problem  and  this  file  requires  7 
conversion.  In  addition  to  this  file,  d^ 
additional  data  files  are  required  per  problem; 
these  files  require  6  or  8  conversion  and  con¬ 
tain  the  following  for  each  detector: 

1.  Detector  coordinates  (in  cm): 

X,  y,  z  If  CT  r  1,  or 

X,  d,  R  If  CT  =  2. 

2.  Program  option  parameters  X  ,  Y  ,  and  Q 
(as  defined  In  3.7.1-D) 

Coordinate  system  parameters  (*10) 

This  data  file  requires  6  or  8  conversion,  and  contains 
the  following  parameters: 

1.  Alphabetic  identification,  SYSTEM,  in  Columns  2 
through  7  of  the  first  card  of  the  file) 

2.  The  following  data  must  be  Input  for  each  of  the 
coordinate  systems  (except  the  base  system,  l.e., 

CS  =  0)  used  to  define  source  points,  detectors. 


or  volumes 


a.  COORN«  coordinate  system  number 

b.  Xq,  coordinates  of  the  origin  of 

the  system  (in  cm)  In  terms  of  the  base 
system. 

c.  X  ,  Y  ,  Z  ,  coordinates  of  a  point  on  the 

X  X  X 

X  axis  of  the  system  10  cm  from  the  origin 
In  terms  of  the  base  system. 

d.  X  ,  Y  ,  Z  ,  coordinates  of  a  point  on  the 

y  y  y 

-fy  axis  of  the  system  10  cm  from  the  origin 
In  terms  of  the  base  system. 

3.8  Card  Identification  Field  for  Problem  Decks 

:.  Columns  63  through  66  are  left  blank.  (These  are  to 
contain  the  Computing  Lab  Job  number  as  specified  by 
the  Computing  Lab.) 

2. .  Columns  6?  and  68  contain  the  problem  deck  number^ 
beginning  with  01  for  any  given  set  of  pz*oblem8. 

3.  Columns  69  through  72  contain  the  card  number^  starting 
with  0001  for  any  given  problem  deck. 

'4.  Column  73  contains  che  last  digit  of  the  year,  l.e., 

0  for  i960,  1  for  1961,  etc. 

3.  Columns  78  through  80  contain  CI7  or  L63. 
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APPENDIX  A 


NOMENCLATURE  AND  SYMBOLISM  FOR  GEOMETRY 
ROUTINE  PLOW  DIAGRAMS 


A-1  Vector  Notation  and  Logic  Symbollam 

Vector  Notation 

1.  All  vectors  are  (3x1)  column  vectors  except  A^^,  A  , 

y 

and  A  which  are  (1x3)  row  vectors. 

Z 

2.  All  primed  vectors  are  In  terms  of  the  base  coordinate 
system;  unprlmed  vectors  are  related  In  the  coordinate 
system  of  the  volume  being  tested. 


3.  For  two  vectors:  A  r 

fx  1 

*a 

and  B  = 

*b 

ya 

yb 

s. 

^b. 

=  >l(^a-Xb)^  (ya-yb)^  (*a“^b^^ 


AxB 


^a*b  '  2bX^ 

Vb  -  ^b^a 


Logic  Symbolism 

— >  Is  read  as  "is  replaced  by">»i.e.^,  a->  b  Is  read 
as  "a  Is  replaced  by  b." 

2.  Branching  relationship  is  defined  as  in  the  figure 
below; 


If  a<b,  the  program  goes  to  block  A 
a  r  b,  the  program  goes  to  block  B 
a>b,  the  program  goes  to  block  C 
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A“2  Test  Function  and  Nomenclature  for  the  Slmple*»Qeometry 
Routine 

Test  Function 

The  element  test  function  used  In  this  routine  Is  a 
linear  Interpolation  on  one  variable.  This  function,  for  the 
argument  Rq,  has  the  form 

-h  -X/t/./  (A-1) 

where  x  and  R  and  x  and  R  ,  are  the  x-value  and  outer 
n  o,n,  n-1  o,n-l 

radius  associated  with  the  n^*’  x-plane,  and  the  n-1®^  x-plane, 
respectively. 

The  function  for  arguments 

Is  similar  to  Equation  A-1  with  the  exception  that  the  new  argu¬ 
ment  Is  substituted  for  R  . 

o 

Nomenclature 

B  Is  gamma  differential  energy  spectra  control  number 

C  Is  Index:  C  ^  0,  volume  search:  determination  of  volume 

In  which  P  lies 

C  =  1,  element  search:  determination  of  element 
In  which  P  lies 

C  =  2,  boundary  search:  determination  of  boundary 
Intercept  B'  for  the  volume 

D  Is  source-detector  distance 

d  Is  subscript  -  refers  to  detector  or  detector-side 

boundary  of  a  volume 

cl  la  Index  of  detector  point  In  set 

f  Is  Index:  f  =  0,  stepping  point  Is  on  an  x-plane 

f  =  1,  stepping  point  Is  not  on  an  x-plane 
f  =  2,  stepping  point  Is  In  base  material 
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1  is  index  of  non-void  segments  such  that  Mj  ^  _1 

numbered  so  that  the  closest  non-void  segment  ' 
to  the  source  point  has  i  =  1. 

K  is  geometry  type  number*  ^  r  0,  cylindric 

=  1,  Cartesian 


K  is  stepping-point  increment 
X  is  direction  cosine  vector  for  P 


A#  A  are  components  of  X,  and  P  direction  cosines  for  the 
^  ^  X,  y,  and  z  directions,  respectively. 


M  is  material  number  (M  r  0  is  void) 

^is  reference  material  number  (listed  in  *8  data  file) 

m  is  subscript  for  maximum 

N  is  volume  number  (input  data) 

Nj^  is  number  of  reference  materials  (see  Sec.  3*1*7-A4) 

n  is  index  of  x-planes  of  a  volume  (n  =  1  denotes  x-plane 
with  smallest  x  value) 

P  is  stepping  point 

?  is  coordinates  of  stepping  point 

,R^  are  x-plane  bounds 

r  is  coordinate,  cylindric  geometry  (r  =  y  ^  z  ) 

S  is  source-point  coordinates 

s  is  subscript  -  refers  to  source  point  or  source-side 

boundary  of  a  volume 

is  index  of  source  point  in  set 

s*  is  index  of  source  set  in  group 

Tj^  is  accumulated  source-to-1  -boundary-intercept  distance 

t^  is  void  thickness  between  stepping  point  and  last  com¬ 
pleted  segment 


U  Is  gamma  calculation  option  (Sec.  3»1.7“A5) 
V  Is  reference  material  counter 
W  Is  segment  length 


Y,Z  are  X“plane  bounds  (Pig.  2-1) 

,y,z  are  cartesian  coordinates  and  vector  elements 
^  Is  material  density 
^  Is  element  test  function 

Test  Functions  and  Nomenclature  for  the  ComDlex-Qeometry 


ROJlSlflg. 

Test  Functions 

1.  For  cyllndrlc  and  rectilinear  volumes:  The  element  test 
functions  used  In  this  routine  are  linear  Interpolation 
formulae.  represents  an  Interpolation  on  one  variable, 
and  ^  represents  an  Interpolation  on  two  variables. 
These  formulae  are  shown  below  for  representative  argu¬ 


ments  . 


2.  For  spherical  volumes:  The  element  test  functions  used 
In  this  routine  represent  a  second  order  Interpolation. 
The  functions  for  a  typical  argument  (representing  R^) 
are: 

V 


Nomenclature 


A  is  coordinate  transformation  matrix  (see  Pig.  2-4  for 
elements  of  A) 

B  is  coordinates  of  the  boundary  Intercept  being  found 

5g  is  coordinates  of  the  source-side  boundary  intercept  of 
the  volume  being  considered 

B(jq  is  coordinates  of  the  detector-side  boundary  intercept 
of  the  volume  containing  the  volume  being  considered 

C  is  index:  C  =  volume  search  (determination  of  volume 

in  which  P  lies) 

C  z  1«  element  search  (determination  of  element 
in  which  P  lies) 

C  z  2,  boundary  search  (ue termination  of  boundary 
Intercept  of  the  volume) 

CS  is  coordinate  system  identification  number  (CS  z  0 
denotes  the  base  coordinate  system)  Input  data 

CT  =  0,  Simple  Cartesian 

1,  Complex  Cartesian 

2,  Simple  cvlinder 
3t  Complex  c /Under 

4,  Compley-rectlllnear  cylinder 
5»  Comblnaolon  I  cylinder 
6,  Combination  II  cylinder 

D  is  source -detector  distance 

13  is  coordinates  of  detector  point  (input  data) 

d(as  subscript’)  refers  to  detector  or  detector-side 
boundary  of  a  volume 

d  is  index  of  detector  point  in  set 

d*  is  index  of  detector  set  in  group 

E  is  element 

f  is  index:  f  >  -1^  P  was  in  the  base  material  but  found  a 

point  in  a  defined  volume 
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f  =  0,  P  not  known  to  be  In  the  base  material 

fa  1,  P  is  In  the  base  material,  JL,  /  0, 

being  found 

f  -  2,  p  is  in  the  base  material,  i  -  0,  0(?*) 

being  found 

f  =  J?  is  in  the  base  material,  z  0,  j2f(P**) 

being  found 

f  P  is  in  the  base  material,  ^  0,  0iF) 

being  found 

f  :  5»  P  is  in  the  base  material,  jt,  =  0,  either 

or  jlf* 

( ^  means* does  not  correspond  tc*^) 

i*  =  6,  P  is  in  the  base  material,  A  /  0,  either 

or  ^ 

g  is  index:  g  =  0,  P  is  not  known  to  be  on  an  x-plane 
g  =  1,  P  is  on  an  x-plane 

1  is  index  of  non-void  segments  such  that  numbered 

so  that  the  closest  non-void  segment  of  the's6urce 
point  has  1=1. 

J  is  index:  J  =  0,  no  volume  imbedded  in  the  base  material 

Impinges  upon  the  source -detector  line 

J  =  1,  some  volume  Imbedded  in  the  base  material 
impinges  upon  the  source-detector  line 

K  is  stepping  point  Increment 

L  is  direction  cosine  vector  for  P 

components  of  X,  and  P  direction  cosines  for  the  x,  y, 
^  and  z  directions,  respectively 

M  is  material  number  (M  =  0  is  void) 

m  is  subscript  for  maximum 

N  is  volume  number  (input  data) 
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n  Is  Index  of  x-planes  of  a  volume  (n  =  1  denotes  x-plane 
with  smallest  x  value) 

0  Is  coordinates  of  the  origin  of  a  coordinate  system  In 
terms  of  the  base  coordinate  system  (input  data) 

P  Is  stepping  point 

P  Is  coordinates  of  stepping  point 

Is  points  used  In  determining  whether  the  source-detector 
line  passes  through  a  subvolume 

q  Is  volume  type:  q  =  0,  master  region 

q  =  1,  region 

q  -  2,  subregion 

Is  element  test  function  (spherical  option) 

R^n'^o  x-plane  bounds  (Pig.  2-1) 

r  Is  coordinate^  cyllndrlc  geometry  (r^  =  y^-f-  z^) 

S  Is  source-point  coordinates 

s  (as  subscript)  refers  to  source  point  or  source-side 

boundary  of  a  volume 

s  Is  Index  of  source  point  In  set 

s*  Is  Index  of  source  set  in  group 

T  Is  accumulated  source- to-boundary  Intercept  distance 

t  is  material  thickness  (along  source-detector  line) 

u  Is  test  parameter:  u  z  y^»  CT<2 

u  :  rp,  CT>1 

ui«U2  are  element  test  parameters 
U^*4U2*  are  envelope  test  parameters 

V  Is  test  parameter;  v  z  CT  <  2 

V  ;  tan“^  ^  ^ 
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are  element  test  parameters 
V^*,V2*  are  envelope  test  parameters 
W  Is  segment  length 

X*  is  element  test  function  (spherical  option) 

Y,Z  are  x-plane  bounds  (Pig.  2-1) 
x,y«z  are  Cartesian  coordinates,  and  vector  elements 
are  x-plane  numbers 
Is  x-plane  test  function 

^  Is  contalnlng-volume  test  function;  Is  the  portion  of 

source-detector  line  bounded  by  the  source  point  and 
the  detector-side  boundary  Intercept  of  the  contain¬ 
ing  volume  “D ) 

^  Is  Index:  ^  z  0,  element  test  function  being  computed 

tor)^ 

^  z  1*  element  test  function  being  computed 

forx^ 

0  Is  coordinate,  cyllndrlc  geometry:  0  z  tan"^  z/y 
Is  Index:  yj  z  0, 

q 

a)  z  1,  IB-S)> 

ZZ  Is  element  test  function 
0  Is  x-plane  bounds  (Fig.  2-1) 

P  Is  material  density 
^  Is  Index: 

Cyllndrlc  Geometry  Cartesian  Geometry 

Working  on  Y  or  r: 

s  0,  mln(R^n)  5  rp  jc  max(RQ)  mln(Y)  f  Yp^  max(Y) 

A  =1,  rp<  mln(R^^)  yp<mln(Y) 

^2=  2,  rp>  max(R^)  yp>  max(Y) 
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Working  on  Z  or  ©: 


mln(  ©p  ^  maxjS^g) 

mln(Z)  $  Zp  < 

max(Z) 

©p  ^  mln(R^j^) 

mln(Z) 

©p  >  maxCR^) 

max(Z) 

APPENDIX  B 


NEUTRON  PLUX-TO-HEAT  CONVERSION  FACTORS 

The  neutron  flux-to-heat  conversion  factors  used  for  these 
programs  must  be  the  heat  generated  per  Incident  neutron  In  a 
material  of  unit  density,  and,  hence,  these  factors  will  have  units 
of  Joules-cm^/gm. 

The  heat  source  due  to  neutron  reactions  with  nuclei  of  the 
target  material  may  be  grouped  according  to  the  manner  In  which 
heat  Is  assumed  to  be  generated  by  the  various  reactions.  For  these 
programs  It  Is  necessary  to  consider  two  modes  of  heat  generation, 
namely,  heat  generated  by  secondary  (or  Induced)  gamma  radiation 
and  heat  generated  by  all  other  reactions  or  portions  of  reactions. 

(As  an  example  of  this  grouping,  consider  the  Inelastic  scattering 
of  neutrons.  It  Is  assumed  that  the  recoil  energy  of  the  scattering 
nucleus  Is  dissipated  as  heat  by  the  second  mode  and  the  de-excltatlon 
photons  generate  heat  by  the  first  mode.) 

The  heating  rates  due  to  secondary-gamma  radiation  cannot  be 
computed  directly.  Using  these  programs,  rather.  It  Is  necessary 
first  to  determine  the  neutron  flux  at  a  representative  number  of 
points  In  the  system  and,  considering  all  possible  sources  of 
secondary  photons,  convert  these  neutrons  Into  secondary- gamma-ray 
sources  using  the  methods  of  Reference  1.  One  may  then  use  these 
sources  together  with  one  of  these  programs  to  compute  the  heat 
generation  rate  due  to  this  phenomenon. 
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The  heat  generation  rate  In  a  given  material  due  to  the  second 
mode  may  be  caused  by  any  one  or  all  of  the  following  reactions: 

•  elastic  scattering, 

•  charged-particle  reactions  (usually  (n,p)  or  (n,d), 

•  particle  emission  from  the  decay  of  activated  residual  nuclei 
( usually  /^*or ^') , 

•  Inelastic  scattering,  and 

•  (n,2n)  reaction. 

The  flux-to-heat  conversion  factor  per  Incident  neutron,  as 
required  for  the  programs.  Is  given  by 

where 

£  Is  the  group-averaged  neutron  cross  section  for  the  reaction 
(in  cm  /gm),  and 

E  Is  the  average  energy  dissipated  as  heat  per  reaction 
(in  Joules). 

For  the  charged-particle  reactions  (from  Ref.  9)t 

I  =  E„+4  -£ 

where  la  the  Initial  neutron  energy, 

Q  Is  the  Q  value  of  the  reaction,  and 
^  Is  the  energy  given  off  as  secondary  photons. 

The  heat  generation  rates  due  to  elastic  and  Inelastic 
scattering  and  the  (n,2n)  reaction  are  derived  In  the  following 


sections 


B-1  Elastic  Scattering 

The  average  energy  loss  per  elastic  scatter,  E,  as  derived 
In  Reference  9,  Is  for  Isotropic  scattering  In  the  center-of-mass 
coordinate  system.  This  assumption  of  Isotropic  elastic  scattering 
Is  not  valid  for  high  neutron  energies  or  for  materials  of  medlum- 
to  high  mass  numbers.  Hence,  a  more  exact  determination  of  the 
average  energy  loss  E  Is  required  In  order  to  predict  heating 
rates  due  to  fast  neutrons. 

The  derivation  of  1  Is  given  In  Section  1.2.1  of  Reference  95 
from  this  derivation,  E  Is  given  by 


^  X'  tPCi)ct£ 
~t£. 


X  '  'PCO^^ 


(B-1) 


where  E^  Is  the  Initial  neutron  energy, 

^  z  2A/(A  +  1)^  (where  A  =  atomic  weight),  and 
p(E)  Is  the  energy-loss  distribution  of  the  degraded  neutrons. 
The  energy  of  the  nucleus  after  collision  Is  related  to  the 
splattering  angle,  0,  by 

^  (B-2) 


Hence,  the  angular  distribution  of  the  scattered  neutrons  may  be 
used  In  the  moments  equation  to  replace  p(E);  or,  since  E  Is  the 
quotient  of  two  moments,  p(E)  may  be  related  to  the  angular  distri¬ 
bution  of  the  neutron  scattering  cross  section.  Certain  of  these 
distributions  for  various  neutron  energies  and  target  materials 
are  given  in  Reference  15. 
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The  equation  for  E  (Eq.  B-l)  may  be  rewritten  in  terms  of 
the  scattering  angle  S  and  the  elastic  scattering  cross  section 
^(Ei,©)  as-  ^ 

Jq  c/-  m  o)  r-C ej ^  cc^^:> 

r-  C£,,  ©)  e) 


>l»  O' )  a£ 

H- 


-7/£, 


or 


F  - 17-  i  ^  -U  i.  ~l  p  £ 

L  -I  ^ 

where  Z  =  cos  Q . 

The  cross  sections  of  Reference  6  were  plotted  against 
Z  and  Equation  B-3  was  used  to  evaluate  E  from  this  data.  For 
the  elements  of  lithium^  zirconium,  tantalum  and  bismuth,  the 
Integral  was  evaluated  numerically  using  Legendre -Gaussian 
quadrature.  The  results  are  shown  In  Figure  B-l. 

Reference  16  Is  a  compendium  of  cross  sections  for  Iron, 
silicon,  aluminum,  and  oxygen.  Among  other  things,  this  report 
lists  coefficients  for  Legendre  expansions  of  the  scattering 
cross  sections;  these  expansions  az^e  of  the  form 


rcr,,e)  ‘  r  czL*nSu 

with^^  defined  by 

X  '-C£„ei Cc^e:>^Cc^o). 

By  definition,  so  that 

<rce,:>  =  2  rrx'  r  c£’„  , 

and 

ir  c£„  e;  eyoL  cc^  s> 

r  cE,,e)  ^cc^e:) 


13^ 


NK:  I4,S3I 


E(M«y) 

PIOURI  l-l.  AVIRAOI  NIUTRON  INIRGY  UMt  POft 
VARIOUS  Z  NUMRIRS 
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This  is  the  integral  term  of  Equation  B-3  where 

Values  of ^  for  iron,  aluminum,  and  oxygen,  computed  from 
data  in  Reference  16,  are  shown  in  Figures  B-2  through  B-4. 

The  computed  values  for -/-show  that  the  actual  heating 
rates  due  to  elastic  scattering  may  be  down  by  a  factor  of  10 
from  those  computed  assuming  isotropic  scattering,  and,  further, 
the  effect  is  most  noticeable  at  high  energy  and  Z  number,  as 
expected  from  the  cross  sections. 

B-2  Inelastic  Scattering  and  the  (n.2n)  Reaction 

The  kinetic  energy  transfer  to  the  residual  nucleus  during 
an  (n,n*)  or  (n,2n)  reaction  is  derived  from  energy  and  momentum 
conservation.  The  final  expressions  are  in  terms  of  the  initial 
neutron  energy  and  either  the  excitation  energy  of  the  residual 
nucleus  or  the  degraded  neutron  energies.  The  expressions  in 
terms  of  initial  neutron  energy  and  excitation  energy  should  be 
used  in  all  possible  cases,  l.e.,  whenever  the  energies  of  the 
resulting  photons  are  known.  If  the  excitation  energy  is  not 
known,  it  must  be  approximated,  and  this  is  done  by  approximating 
the  energies  of  degraded  neutrons  (Sec.  B-2. 3)  and  finding  the 
recoil  energy  from  an  application  of  the  conservation  of  energy. 

Section  B-2. 4  shows  that  for  the  energy  range  of  Interest, 
the  excitation  energy  (and  also  mass)  is  essentially  the  same  in 
both  the  laboratory  and  center-of-mass  coordinate  systems. 
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B-2 , 1  The  (n^n*)  Reaction 


Consider  the  neutron-nucleus  system  shown  in  the  sketch. 


L  system: 
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O^- 
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before 


CM  system: 
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after 


Prom  conservation  of  energy  and  momentum  In  the  CM  system, 

Va  ,  Ci/;-  '/2  /ii/^  =  */2 

(B-4) 

where  m  Is  the  mass  of  the  neutron, 

M  Is  the  mass  of  the  nucleus, 

€.  Is  the  excitation  energy  left  with  the  nucleus,  and 
the  i>r's  are  as  shown  In  the  sketch. 

The  momentum  of  the  neutron  before  the  reaction,  as  measured 
In  the  L  system,  must  be  equal  to  the  momentum  of  the  CM  system  as 
measured  In  the  L  system;  hence, 

//O/T  =  QM  • 
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The  velocities  V’2  and  are  related  by  the  cosine  law: 


Km  ^ 


(B-6) 


This  relationship  is  shown  In  the  sketch  below,  which  Is  a 
superposition  of  significant  portions  of  the  L  and  CM  diagrams 
sketched  above. 

"^2 


Solving  Equations  B-4,  B-3t  and  B-6  for  vl  in  terms  of 

gives 

L .A  . 

for  Eq  2  i  />7t  cT;^, 

A2  •i^  ,  and 
in 

'  (A  +  1)2 

Energy  Is  conserved  In  the  L  system  and,  therefore. 


E^  =  E'-f  , 


for 


£  £  Vt  ^  lt*  j 

a  /1 ^ 


(B-8) 
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And  from  Equations  B-7  and  B-8, 


Prom  Equation  B-9-i 


ote  ,  . 

CxiO  & 

1 

Also,  for  constant  and  ^ 

~  ^0  /  C/-  )  ~yt*/ 

--  2  / 

The  average  value  of  Ej^  is  taken  to  be 

/^Myiro  ig^  4 

It  =  yC-A  pee)  ^  <9^4, 

where  p(©)  la  the  angular  scattering  probability  function  or, 
equivalently  (for  this  purpose),  the  angular  dependent  cross 
section. 


If  isotropic  scattering  In  the  center-of-mass  system  Is 
assumed. 


* 


4; 


S, 


U£ 


or 


^  -  E  i — ^ 

^of  } 

or,  in  terms  of  Eq  and  E, 


i  =  E 
^  ^(MO  ® 


£; 

>9 


(B-lOa) 


(B-lOb) 
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B-2.2  The  (n,2n)  Reaction 

The  equations  describing  the  emission  of  the  first  neutron 
are  similar  to  those  describing  the  (n«n')  reaction* -so  that  It 
may  be  assumed  that  the  kinetic  energy  of  the  nucleus  before 


emission  of  the  second  neutron*  as  measured  In  the  L  system 
(see  sketch)  Is  given  by  Equation  B-10. 
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Prom  the  conservation  of  energy  and  momentum  in  the  CM  system* 

where  M*  Is  the  mass  of  the  residual  nucleus  and*  to  a  good  approxi¬ 
mation*  is  given  by  M*  z  M-m  (except  for  calculating  D)f 
S'iB  the  excitation  enex^  of  the  residual  nucleus; 

D  is  the  mass  defect  and  is  given  by  D  =  (M-M*-m)c  *  and 
all  other  parameters  are  as  defined  above. 
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By  superimposing  portions  of  the  L'  and  CM'  diagrams  sketched 


above,  it  may  be  seen  c/J  ^  are  related  by 


(B-12) 


Solving  for  '/z /m  tr*  in  terms  of 


gives 


(B-13) 


for 


Conservation  of  energy  In  the  L'  system  requires  that 

r)',£-V£'-£> 

£;  =  y*. 


Prom  Equation  (B-l4) 


144 


Also,  for  constant  £j 


j 


,  and 


Analogous  to  for  Isotropic  scattering,  the  average  value  of 

or 


S^-'EdE  _p  .  „ 

'  '■W  7/-  ■  <V  J 


F.  „  ,  ®i£i£: 


(B-15a) 


or.  In  terms  of  E'  and  E",  the  energies  of  the  first  and  second 


neutrons  emitted,  and  E^  the  Initial  neutron  energy. 


ACA*0 


(B-15b) 


B-2.3  Average  Value  of  the  Degraded  Heutron  Energy 

Where  a  large  number  of  closely  spaced  energy  levels  are 
Involved,  the  methods  of  statistical  mechanics  may  be  used  to 
determine  the  value  of  the  degraded  neutron  energy.  Thus,  this 
theory  will  not  represent  the  high  end  of  the  degraded  neutron 
spectrum  very  well  since  this  corresponds  to  leaving  the  nucleus 
In  a  low  excited  state  where  the  levels  are  few  and  far  apart. 
Further,  In  lead.  Iron,  and  most  light  elements,  the  level  spacing 
Is  large  so  as  to  make  this  method  a  rough  approximation  at  best. 
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The  differential  probability  of  emission  of  a  neutron  with 
energy  E'  from  a  compound  nucleus  with  energy  E  is 

where  ©  is  the  nuclear  temperature  and 
k  is  a  normalization  constant. 

The  average  value  of  E*  is  taken  to  be 


where,  for  the  (n,n')  reaction, 

E»*  is  the  initial  neutron  energy  and 
E*  is  zero; 

for  the  (n,2n)  reaction,  and  for  the  emission  of  the 
first  neutron, 

E**  is  the  initial  neutron  energy,  and 

E*  is  the  energy  such  that  the  second  neutron  is  emitted 
with  zero  energy  and  the  nucleus  is  left  in  the  ground 
state  (this  is  the  binding  energy  of  the  remaining 
neutron  in  the  nucleus  after  the  first  neutron  has  been 
emitted); 

for  the  emission  of  the  second  neutron, 

E**  is  the  initial  neutron  energy  less  that  energy  required 
to  emit  the  first  neutron  with  zero  energy  and  leave 
the  nucleus  in  the  ground  state  (this  is  the  initial 
neutron  energy  less  the  binding  energy  of  the  last 
neutron  in  the  compound  nucleus)  and 
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E*  is  zero. 


Therefore, 


/-  m^ti 

1--^ — 


OC  =  £■  .and  ^  /&  ■ 


(B-16) 
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APPENDIX  C 


NEUTRON  REFERENCE  MATERIAL  COMPARISON 

To  assist  In  the  selection  of  base  materials  to  be  used  In 
these  IBM  programs,  a  series  of  problems  have  been  run  using  C-17 
In  a  simple  cylindrical  geometry.  The  reactor  core  was  composed 
largely  of  carbon,  and  various  shielding  materials  were  placed 
adjacent  to  the  end  of  the  core.  For  a  shielding  material  of 
lithium  hydride,  detector  points  were  located  at  various  distances 
from  the  face  of  the  reactor  In  the  lithium  hydride.  Eight 
different  base  materials  were  used  In  computing  the  dose  rate, 
spectra,  and  heating  rate. 

Figure  C-1  shows  a  comparison  of  heating  rates  as  calculated 
for  the  different  reference  materials.  Figure  C-2  shows  a  compari¬ 
son  of  dose  rates  as  calculated  for  the  different  reference  materials. 

In  order  to  get  some  Idea  as  to  the  reliability  of  the  data 
calculated  by  these  IBM  programs  additional  calculations  have  been 
made  of  the  neutron  differential  number  spectra  and  compared  to 
those  calculated  by  a  multigroup  multiregion  diffusion  code. 

These  comparisons  are  shown  In  Figures  C-3»  C-4,  C-5*  C-6,  and  C-7. 

It  Is  apparent  that  for  "O"  distance  from  the  core  face  the  non- 
hydrogenous  reference  materials  differential  number  spectra 
calculations  do  not  agree  with  those  calculated  by  the  diffusion 
code.  This  Is  due  to  the  Inability  to  fit  the  moments  data 
accurately  for  penetration  distances  of  less  than  10  cm^/gm. 
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FIGURE  C-1.  HEAT  DEPOSITION  RATES  IN  LITHIUM  HYDRIDE 
FOR  VARIOUS  BASE  MATERIALS 
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FAST-NEUTRON  DOSE  RATE  (mrem/hr-wa«) 


NPC  l4,S3i 
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FIGURE  C-2.  FAST-NIUTRON  DOSE  RATES  IN  LITHIUM  HYDRIDE 
FOR  VARIOUS  BASE  MAHRIALS 
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NEUTRON  FLUX  (neutrons/cm^-sec-Mev-waff) 


L52 
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riGURI  C-4.  COMPARISON  OF  C-17  AND  ZOOM  SPICTRA 
FOR  SORATfD  CARRON 
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NEUTRON  FLUX  (neutrons/cm^-sec-watt-Mev) 


ENERGY  (Mev) 


FIGURE  C-5.  COMPARISON  OF  C-17  AND  ZOOM  SPECTRA 
FOR  BORATED  BERYLLIUM  OXIDE 
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NEUTRON  FLUX  (neutron$/cm^-sec-Mev>watf) 
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FIGURE  C-6.  COMPARISON  OF  C-17  AND  ZOOM  SPECTRA 
FOR  BORATED  BERYLLIUM 
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NEUTRON  FLUX  (neutrons/cmJ-sec-Mev-watt) 


ENERGY  (Mev) 


FIGURE  C-7.  COMPARISON  OF  C-17  AND  ZOOM  SPECTRA  FOR 
BORATED  POLYETHYLENE  PLUS  FIBER  GLASS 
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APPENDIX  D 


PROGRAM  LIBRARY  DATA 

The  Library-Types  1  and  2  currently  being  used  with  programs 
C--i7  and  L-63  are  presented  in  this  appendix.  The  data  are  listed 
directly  from  the  data  cards j  however,  captions  have  been  added 
to  each  Library-Type  1  material  deck.  The  formats  for  the  Library- 
Type  1  decks  and  the  Library-Type  2  deck  are  shown  in  Figures  3-1 
(Sec.  3*2)  and  3-2  (Sec.  3*3)»  respectively. 
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LIBRARY-TYPE  1 


MATERIAL  LITHIUM  HYDRIDE 


NEUTRON  ENEROjY 
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8  -335684-3+1 251  07+0-684 140-4+490505-1 -559075-5-5 86205+1 
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